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ABSTRACT: We study the form factors of the Konishi operator, the prime example of non- 
protected operators in M = 4 SYM theory, via the on-shell unitarity method. Since the 
Konishi operator is not protected by supersymmetry, its form factors share many features 
with amplitudes in QCD, such as the occurrence of rational terms and of UV divergences 
that require renormalization. A subtle point is that this operator depends on the spacetime 
dimension. This requires a modification when calculating its form factors via the on-shell 
unitarity method. We derive a rigorous prescription that implements this modification to all 
loop orders and obtain the two-point form factor up to two-loop order and the three-point 
form factor to one-loop order. From these form factors, we construct an IR-finite cross- 
section-type quantity, namely the inclusive decay rate of the (off-shell) Konishi operator 
to any final (on-shell) state. Via the optical theorem, it is connected to the imaginary part 
of the two-point correlation function. We extract the Konishi anomalous dimension up to 
two-loop order from it. 
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1 Introduction 


So far, the framework of quantum field theories (QFTs) is very successful in describing 
the high-energy processes measured at colliders such as the LHC. However, theoretical 
predictions are usually restricted to the weak-coupling regime, which admits a perturbative 
expansion in terms of the small coupling constants. The individual contributions to the 
perturbation series can be calculated via Feynman diagrams. Thereby, a large proliferation 
of diagrams is in general encountered when one proceeds to higher-order corrections, and 
hence concrete calculations are mainly restricted to the first few orders. 

The investigation of alternative techniques that bypass this limitation is thus of high 
importance. It might not only allow to push perturbation theory to higher orders, but could 
also deepen our understanding of the fundamental principles and mechanisms encoded in 
QFTs. The so-called ‘on-shell’ techniques are such an alternative. They allow one to 
build amplitudes from simpler amplitudes with a lower number of external legs and loops 
via recursion relations |1, 2] and unitarity [3, 4]. They have been successfully used in 
supersymmetric gauge theories as well as in QCD, see [5-7] for pedagogical reviews and 
references therein. 

In particular, the maximally supersymmetric Yang-Mills (V = 4 SYM) theory with 
gauge group SU(N.) in four dimensions has played an important role in the aforementioned 
developments. According to the AdS/CFT correspondence [8-10], it has a dual description 
in terms of a string theory, allowing its study also at strong coupling. Moreover, in the 
planar limit [11], it shows signs of integrability at weak as well as at strong coupling, which 
is believed to be present even at any coupling. Based on the conjectured integrability, new 
predictions for the spectrum, i.e. for the anomalous scaling dimensions of gauge-invariant 
composite operators, were made; see [12] for a review. This rises the hope that the theory 
is exactly solvable, and it is hence sometimes even referred to as the “harmonic oscillator 
of the 21% century”. 

Given the success of the aforementioned on-shell techniques for amplitudes, it is an 
intriguing question whether they can be applied for determining off-shell quantities such 
as correlation functions or the anomalous dimensions as well. A bridge between the purely 
on-shell amplitudes and the purely off-shell correlation functions is provided by form fac- 
tors. In particular, they also contain the information necessary to determine the anomalous 
dimensions. An n-point form factor describes the overlap of an off-shell initial state, de- 
scribed by a composite operator, into an on-shell final state consisting of n elementary 
fields. It is given by 


Fo(1,.-.,n) = f dP ae #*(1--- nJO(a)|0) = (21) 8P (q— Y pi) (1 nO 00) , 


i=1 
(1.1) 
where the particles labeled by i = 1,..., carry individual on-shell momenta p; and the 
operator O carries off-shell momentum q. If the number n of the external fields exactly 
matches the number of fields contained in O, the form factor is called minimal. Minimal 
form factors with n = 2 points are denoted as Sudakov form factors. 


In M = 4 SYM theory, the most intensively studied form factors are the ones of the 
half-BPS operator 


Opsps = tr(¢u¢s)) , (1.2) 


where the parentheses denote traceless-symmetrization of the indices J, J = 1,...,Ng of 
the Ng scalar field flavors. This operator belongs to the stress-tensor supermultiplet. Its 
Sudakov form factor was first studied by van Neerven [13] and analyzed up to four loops 
[14, 15] in the recent past. The Sudakov form factor exhibits exponentiation [16-18], a 
feature which was seen to be the key for predicting the all-loop IR behavior of scattering 
amplitudes [19]. 

The form factors of the stress-tensor multiplet with general n external legs can be 
analyzed in analogy to scattering amplitudes with modern on-shell techniques. The n- 
point form factor of the bosonic operator (1.2) was first studied in [20, 21], and later 
generalized to the full stress-tensor multiplet in [22, 23]. Up to one-loop order, compact 
expressions for general n-point maximally-helicity-violating (MHV) as well as some next-to- 
MHV (NMHV) form factors have been computed in [20, 22-26]. The two-loop three-point 
form factor was computed in [27]. The form factors of half-BPS operators with k scalar 
fields, as well as the corresponding supermultiplets, have been studied in [21, 28, 29]; n- 
point tree and one-loop MHV results are presented in [28] and the mininal form factors 
(for n = k) were computed at two-loop [29]. Form factors have also been studied at strong 
coupling via the AdS/CFT correspondence [30], and a Y-system formulation was given in 
[31] for AdS3 and in [32] for AdS;. 

The aforementioned studies have shown that form factors share very similar recursive 
and analytic properties with scattering amplitudes, at least for the protected operators. 
Moreover, the robust set of on-shell techniques for computing on-shell objects is also ap- 
plicable in this case. This rises the hope that also fully off-shell quantities can be studied 
using on-shell methods, and that such an enhancement of the toolkit allows to detect new 
features of the theory. Indeed, it was found that certain correlation functions can be con- 
structed via generalized unitarity from amplitudes, form factors and their generalizations 
involving several operator insertions [25]. In the recent parallel work [33], one of us has de- 
termined at tree level the minimal form factors of a generic operator and at one-loop order 
their cut-constructible parts. The one-loop results yield the complete one-loop dilatation 
operator of the theory. 

Scattering amplitudes as well as form factors are themselves not physical observables, 
since they contain infrared (IR) divergences from the integration of loop momenta. Adding 
the so-called bremsstrahlung contributions, their IR divergences from the real emissions of 
soft and collinear particles cancel the IR divergences coming from virtual loop corrections 
according to the Kinoshita-Lee-Nauenberg theorem [34, 35], and one obtains an observable. 
In particular, the cross sections are free of IR divergences and hence physical observables. 
They are, however, in general not well defined in a CFT such as V = 4 SYM theory, where 
asymptotic states are ill defined. Some cross-section-type quantities have been defined by 
using coherent states as asymptotic states [36]. Alternatively, we can consider the decay 
of an initial off-shell state created by an operator O(q) with timelike momentum (q? > 0) 


into any final on-shell multi-particle state. The probability of this inclusive decay is the 
total decay rate of O(q). This decay process may occur as part of a total cross section of 
a scattering process in which O(q) is produced as an intermediate state.! The probability 
for the inclusive decay of O(q) into a final state X with total momentum q = px is defined 
by 


oola) = >, 8) (q— px) |(X|O(0)|0)/ , (1.3) 
X 


where the sum ensures that the quantity is inclusive, i.e. all contributions, which are 
specified by the number and type of the particles in the final states, are integrated over 
the respective phase space and are summed up. This cross-section-type quantity depends 
on the matrix element (X|O(0)|0), which is precisely the form factor of O with final state 
X. Via the optical theorem, (1.3) is related to the imaginary part of the (time-ordered) 
two-point correlation function (0/O(x)O(0)|0) after transforming to momentum space. 


Finally, although not considered in this paper, we would like to mention that by 
modifying (1.3), ‘event shapes’ such as energy or charge correlation functions were studied 
in NV = 4 SYM theory [25, 38-40]. Also, Wilson coefficients for deep inelastic scattering 
were considered [41]. For simplicity, we will follow the terminology of [39] and denote the 
cross-section-type quantity defined in (1.3) as total cross section, or simply cross section. 


In this paper, we will study the form factor (1.1) and the cross section (1.3) for the 
Konishi operator as a first example for an operator that is not protected by supersymme- 
try. Hence, UV divergences appear in addition to the aforementioned IR divergences that 
already emerge for protected operators. The Konishi primary operator is given by 


K = 6" tr(¢rdy) , (1.4) 


where sums over all J, J = 1,..., Nọ scalar field flavors are implicitly understood. In strictly 


D = 4 dimensions, we have Ng = 6. The Konishi scaling dimension Ax = AO) tn consists 


of the bare dimension A®) = 2 and an anomalous dimension yx. It is a power series in the 


coupling constant 


2 
N, 
g? = PM e (ame) , (1.5) 


which depends on the Yang-Mills coupling constant gym as well as the number of colors Ne 
and is the loop-counting parameter in the modified dimensional reduction (DR) scheme in 


'The operator may be of different physical origin. For example, it can be part of a vertex that couples 
to a massive particle in an effective Lagrangian. Then, (1.3) yields the decay rate of this particle. A 
concrete example from the Standard Model is an effective Higgs-gluon vertex H tr(Fy,F"”) obtained by 
integrating out a heavy quark loop, see e.g. [37]. The operator may also be a (conserved) current describing 
a two-particle scattering. Examples of this type are ee~ annihilation into a virtual photon or Drell-Yan 
scattering, where the two incoming particles are annihilated into a virtual photon or gluon, respectively, 
exciting the QCD vacuum and decaying into quarks, gluons etc. 


D = 4— 2e dimensions.” In the planar limit, the Konishi anomalous dimension is given by? 


yg = 6[29? — 894 + 56g — 16(26 — 663 + 15¢5)g° 


1.6 
+ 16(158 + 72¢3 — 54¢? — 90¢5 + 315¢7)g"9] + O(g!) , oo 


where the one- and two-loop contributions, which we reproduce as a check in this paper, 
were obtained by explicit Feynman diagram calculations in [45, 46] and [47—49].* 

The operator (1.4) is the primary operator of the Konishi supermultiplet. Its anoma- 
lous dimension given in (1.6) was mainly obtained by considering certain descendent op- 
erators within the Konishi multiplet rather than the Konishi primary operator (1.4). This 
is possible since all members of a supermultiplet have the same anomalous dimension.” In 
fact, we will see that the Konishi primary defined in (1.4) and involving a sum over the Ng 
scalar field flavors depends on the dimension D, since Ng = 10 — D is required to ensure 
supersymmetry. This becomes important when regulating the divergences by continuing 
the theory from D = 4 to D = 4 — 2e dimensions. 

We will apply four-dimensional unitarity in order to compute the form factors. Within 
this framework, all on-shell component fields can be conveniently combined into Nair’s 
N = 4 on-shell superfield [65]. The on-shell superfield reads 


A,B A,B,C 


&(p,n) = 94(p)-+n4 walp)+2 x PAB (p) T jo (p)+n'n?n?ntg-(p) , (1-7) 


where ņ^ are Grassmann variables that encode the flavor and helicity of the component 
fields. Pairs of upper and lower SU (4) R-symmetry indices A, B,C, D = 1,...,4 are always 
understood to be summed. In the above superfield, the six real on-shell scalars @; trans- 
forming in the fundamental representation of SO(6) are represented via the anti-symmetric 
product representation of two fundamental SU(4) representations, 64g = ¢1(07) AB, em- 
ploying the isomorphism of the Lie-algebras so(6) and su(4) induced by the o-matrices 
(or)aB = —(01)Ba- 

Using (1.7), each n-point scattering amplitude with fixed total helicity can be efficiently 
packed into a single superamplitude. In analogy, also the form factors for the BPS operator 


?The DR scheme employs dimensional reduction of ten-dimensional M = 1 SYM theory to D = 4 — 2e 
as regularization [42, 43] and for the subtraction of the divergences a modified minimal subtraction which 
absorbs the same finite terms in addition to the UV divergences into the renormalization constant as the 
famous MS scheme [44], leading to the factor (4me~7£)°*. 

3Note that there are no non-planar corrections to yg at the first three loop orders. 

“The Konishi anomalous dimension yx is currently known up to five loops from field theory calculations 
and up to nine loops from the conjectured integrability. The three-loop result was conjectured in [50] and 
confirmed in [51, 52]. The four-loop result was determined by calculating the wrapping corrections to the 
integrability-based asymptotic dilatation operator in [53, 54] and by a computer-based direct calculation 
in [55]. The integrability-based four-loop expression of [56] matches this result. The five-loop result was 
predicted from integrability in [57-59], and confirmed in [60] from an OPE analysis of the four-point 
correlation function of stress-tensor multiplets. The results at six [61], seven [62], eight [63] and nine loops 
[64] are so far only based on the conjectured integrability. 

°Working with certain descendants which are non-singlet states of the SU(4) R-symmetry instead of the 
primary operator (1.4), which is an SU(4) singlet, simplifies the calculations in both, the field theory and 
integrability-based approach. 


(1.2) can be packed into super form factors if the BPS operator is expressed in terms of 
the scalar fields @4p as 


Ozps = tr(¢aBdcp) — Z easco tr(¢ F bpp) , (1.8) 


where the last term subtracts the trace in the space of scalar flavors.°® 


Without loss of generality, we will focus in the rest of this paper on its particular 
component 


Opps = tr(PaBoas) , (1.9) 


where doubled indices are not summed. Expressing also the Konishi operator in terms of 
the scalar fields ¢4p yields 


Ke = enon tr(¢aBpop) = tr(di2¢34) — tr(ġ13%24) + tr(p14%23) , (1.10) 


where the subscript 6 reminds us that the operator is identical to the Konishi primary (1.4) 
only for Ng = 6, i.e. only in strictly D = 4 dimensions. 

There is a subtlety originating from the fact that in D # 4 dimensions the Konishi 
operator K in (1.4) cannot be identified with Xe in (1.10). The four-dimensional unitarity 
method directly applies to the operator Kg. In this formulation, the operator stays the 
same if the encountered IR and UV divergences are regularized by changing the spacetime 
dimension from D = 4 to D = 4 — 2e. But in D = 4 — 2e dimensions the Konishi operator 
K is not identical to the operator Kg. Hence, the unitarity-based results for Kg do not 
directly yield those for the Konishi operator K. Instead, modifications have to be made 
which take into account that one should have used K and not Ke in order to obtain the 
results for the Konishi operator regularized in D = 4 — 2e dimensions. 

In the main part of the paper, we elaborate on the ideas mentioned above. In section 
2, we discuss two-point correlation functions of gauge-invariant local operators, their renor- 
malization and the transformation to momentum space. We identify the imaginary part of 
such a correlation function with the cross section defined in (1.3). Finally, we present the 
general strategy of computing the total cross section for a given operator using its form 
factors as building blocks. 

In section 3, we present our computation of the form factors for Kg at the one- and 
two-loop orders, which are based on the unitarity method and on-shell superspace. Since 
the Konishi operator is not protected, several interesting features appear in the results 
which have not occurred for amplitudes or BPS form factors in M = 4 SYM theory, e.g. 
UV divergences and rational terms. 

In section 4, we discuss in detail the aforementioned subtleties arising from the fact 
that in D = 4 — 2e dimensions the Konishi operator K cannot be identified with Kg. We 
derive a rigorous prescription of how to implement the substitution of Kg by K in the 
results of the previous section and give final results for K. 

In section 5, we present the computation of the cross section starting with the BPS 
operator up to one-loop order as a simple example to make the reader familiar with our 


Note that (01)4°(o7)cB + (03)4° (o1)cB = —26r56“ B, where (o7)4? = $e4?°P (ar)op. 


strategy. We find the expected non-trivial cancelation of the IR divergences between real 
and virtual channels. Then, we compute the cross section for the Konishi operator up to 
two loops. We extract the renormalization constant and hence the anomalous dimension 
from the UV divergence of the bare result. They match the known expressions. We 
present the finite result for the renormalized cross section and discuss its dependence on 
the renormalization scheme. 

Finally, in section 6 we summarize the main results of our paper and the interesting 
features associated with them. We also present some future directions and open questions. 

In the appendices A, B and C, we give some further conventions and explicit results 
for the occurring loop integrals as well as Passarino-Veltman (PV) reduction formulae. 
Appendix D provides some cross checks for the one-loop three-point Konishi form factors. 
In appendix E, we present some details on the phase space integrals occurring in section 5. 
A way to extract the anomalous dimension directly from the two-point Konishi form factor 
is given in appendix F. In appendix G, we discuss the renormalization-scheme dependence 
of the cross section. In the final appendix H, we summarize direct Feynman-diagrammatic 
calculations of the one- and two-loop form factors for the BPS and the Konishi operator, 
which serve as checks for our approach and guided us to the modifications discussed in 
section 4. 


2 Cross sections for two-point correlation functions in a nutshell 


In this section, we review some facts about the form of the two-point correlation function of 
a renormalized composite operator in spacetime and in momentum space. Via the optical 
theorem, its imaginary part yields a cross-section-type quantity. We present our strategy 
of computing this quantity from the form factors of the respective operator. 


2.1 Renormalization of composite operators and their two-point functions 


Gauge-invariant local composite operators can be regarded as external states of M = 4 
SYM theory, and they can occur in correlation functions in the same way as the elementary 
fields. Such correlation functions in general contain UV divergences which are associated 
with the presence of these operators, requiring their renormalization in analogy to that of 
the elementary fields and vertices of the theory. In this paper, we only consider composite 
operators that are eigenstates under renormalization. Such a renormalized operator is 
given in terms of the bare operator as 


OR = Zo(g,6€)Op . (2.1) 


The renormalization constant Zo depends on the coupling constant g and absorbs the UV 
divergences, which appear as poles in e when the theory is regularized by changing the 
spacetime dimension from D = 4 to D = 4— 2e. The renormalization constant determines 
the anomalous dimension 


x D 0 
ro = Yo 9” = lim «95, 108 Zo ; (2.2) 
l=1 


O(0) C> A(x) q G O- 


Figure 1: The two-point function in position space and momentum space. 


which is added to the bare scaling dimension A) in order to obtain the conformal dimen- 


sion Ao. Since yo is finite when the limit € > 0 is taken in the above equation, the form 
of Zo as a power series in g is fixed to 


© 28 (1) (1))2 (2) 
Y 
Zo = exp ( TW) =1 +P + g (Se t + Jo ) + O(g°) . (2.3) 


Conformal symmetry also completely fixes the form of the two-point function of the 
operator Og. In Minkowski spacetime, it reads 


M 


ee — a (0) 
(—a? + i0)40 p270 , Ao = Ao + Yo , (2.4) 


G20,R(£) = (0|OR(2)OR(0)|0) = 
where our conventions for the i0 description are given in appendix A. The parameter u has 
the dimension of mass and is introduced in order to fix the mass dimension of Ggo,R to 


2A), The coupling-dependent dimensionless factor M has a perturbative expansion as 


M= y un, (2.5) 
L=0 


and it can be absorbed into the normalization of Op. 

We will work in momentum space, and hence need the Fourier transformation of (2.4). 
According to appendix A, it is given by 

D 
Grok (7?) = fare et? Goo R(T) = (—i)2P-24072 HG Ao) =) 2 
(2.6) 

When expanding the above expression first for small g and then for small €, one obtains 
4.-poles for any k > 1, which for k > 2 are proportional to powers of yo [66]. Since 
G2o0,R(«) is the finite (renormalized) Green function, these poles cannot come from UV 
divergences. In fact, they arise from integrating over the origin x = 0 of spacetime, 
where G20 R(x) is singular. This can be most easily seen for the half-BPS operator Opps 
defined in (1.2). Since this operator is protected, ygps = 0, and all poles of order k > 2 
disappear, but a simple +-pole remains. In momentum space, this pole is associated with 
the one-loop bubble integral. It is obtained when inserting Fourier expressions for the 
two scalar propagators’ aan connecting the two operators as depicted in figure 1 


"In D dimensions, the scaling dimension of a scalar field is given by AY = 2 —-l=l1-e. 


and performing the integration over x in (2.6), which yields a 6-function of momentum 
conservation. For the tree-level two-point function, the steps are as follows: 


1 aot 1 1 
— r > Eem oo (2.7) 
(—x2 + i0) FT (27)? (l — q)? (~? — i0)€e 
This simple pole (for the BPS operator) and all the further 4.-poles, k > 2, (for non- 
protected operators) are absent when taking the imaginary part of the momentum-space 
Green function (2.6). 

As we will see in the next subsection, via the optical theorem the imaginary part 
of (2.6) yields a cross-section-type quantity: the probability of the inclusive decay of the 
renormalized operator Org. It has to be finite in the limit € — 0, since it is free of IR 
divergences and — due to renormalization — also of UV divergences. 


2.2 Two-point correlation functions and cross sections 


Via the optical theorem, the imaginary part of a two-point correlation function is related 
to the inclusive decay width of the renormalized operator Or with off-shell momentum q, 
where q? > 0. As motivated in the introduction, we will simply denote this as cross section 
co in this paper. It is given by® 


oor = Im[2i Gr0,x(¢)] = J. (27)? 8P) (q — px) |(X|Or (0)|0)/? , (2.8) 
X 
where one sums over all final on-shell states X, and the squared matrix element is given 
by the product of two form factors? 


Fo x = (X|O(0)|0) . (2.9) 


The form factor has the perturbative expansion 
OO 
fox =y g" foy (2.10) 
é=0 


where g is the parameter of the loop expansion. Concretely, in M = 4 SYM theory in 
the modified dimensional reduction (DR) scheme, the coupling constant is given in (1.5). 
Moreover, the summation over all final states X in (2.8) involves in particular a summation 


) 


over the number n of particles in the final state, i.e. of the n-point form factors FY n OVET 


n. The number n is directly related to powers of the Yang-Mills coupling constant gym. In 
analogy to amplitudes (see e.g. [67]), the n-point form factors possess a decomposition in 


terms of the possible color structures as 


(Ee n— a a £ 
FY) Canpa mY) =g, eee eT 2) FS, Apa Notiy}) 
Sn /Zn (2.11) 


+ multi-trace terms , 


SNote that the factor of (—i) appearing in (2.6) due to the Wick rotation (see appendix A) must be 
removed before taking the imaginary part. Hence, we have to take the imaginary part of iGo. 

°Here and in the following, we understand that prefactors (27)?6)(q — px) ensuring momentum 
conservation have been stripped off from the form factors. 


where T°, a = 1,..., N2 — 1, are the gauge-group generators of SU(N.) normalized as 
er) = 0” . (2.12) 


In (2.11), the it? particle, i = 1,...,n, with momentum p; carries the adjoint gauge-group 
index a;. Via Nair’s superfield (1.7), its flavor and helicity are encoded in terms of the 
Grassmann variables ņ;, on which the color-ordered super form factors Pe on the rhs. 
also depend. 

The imaginary part of (2.6) can be obtained by taking the discontinuity, which for 


timelike (q? > 0) momentum reads!” 


271 


r&r — x) 


2i Im (—q* — i0)® = (—q? — 10)” — (—q? + i0)” = (q. (2.13) 


Using this relation in order to determine the imaginary part of (2.6) and then inserting 
the result into (2.8) yields 


cor _ M(g) TAPE - = DO — gr + $45) / g )” (2.14) 
D l 
2 


l) MO NAZ - A®) TAg = 


where we have divided cor by its classical part o% = Imf2i C a Indeed, as 
mentioned at the end of the previous subsection, both o% and co, are free of +-poles, 
since the poles are canceled by the extra I'-functions introduced via (2.13). This can also 
directly be seen for the bubble integral in (2.7): its imaginary part is obtained by applying 
a double-cut, which just yields a finite constant. 


By taking the logarithm of (2.14), we can expose the dependence on q? as follows: 
JOR q 
log 70 = yo log A +C+O(e), (2.15) 


where the constant C is scale-independent but depends on yo and the expansion coefficients 
of the normalization factor (2.5) as 


MQ) 
C=9 Go =(= 2%) 


M(©) oe 
M® 1/MQM\? 3-7 
‘ -= a2 4q (2) 6 
tg (san (To) r o a-m) ) +00 ). 


It is also renormalization-scheme-dependent as discussed at the end of section 5. However, 
the log q? term is universal and scheme-independent. The anomalous dimension is given by 
the coefficient of log 5. In this paper, we will verify this structure for the Konishi operator 
up to two loops. 


1QOur conventions for the i0 description are given in appendix A. 


—10-— 


Strategy of computing cross sections 


The cross section is obtained from (2.8) in more detail as follows: 


-gg z {OC SO} e 
colors spins 
e E 


Mn 


This relation holds for both, the bare and the renormalized cross section, if Fn represents 
the bare and the renormalized form factors, respectively. The evaluation of (2.17) requires 
three main steps: (1) determining the form factors F,, (2) taking the absolute square of 
Fn, and (3) performing the n-particle phase-space integrals. More concretely, (2.17) is 
expanded in powers of g as follows: 


£42 


g= Sog” : = 20 (2-n) fas, MEt- n (2.18) 
é=0 


where the squared matrix elements are given by 


n m £ 
1 m—k 
MO = => Ioh X X a _ ut ({a:, pi, nm} F By, — "(fas pem}) y 
ay i=1 


k=0 1=0 
(2.19) 
in which ÊẸ ({ai, pi; m}) is the -loop n-point non-color-ordered super form factor defined 
n (2.11), and F*({a;,p;, m}) is its complex conjugate.!! Moreover, k in N*MHV is called 
MHV degree, which refers to terms in FO with a specific degree in 7. For the BPS 
and Konishi operator considered in this paper, the MHV form factors have degree 4 in 7 
and m = n — 2 is fixed. The squared matrix element involves sums over all numbers n and 
types of external particles as well as their color degrees of freedom. The sum over the types 
of particles is given in terms of integrations over the fermionic variables né, A = 1,2,3,4, 
and a sum over the MHV degree k. 
Given the squared matrix elements, as a next step, the integration over the phase space 


of the n particles in the final state has to be performed. The respective measure is given 
by 


dPs, = Gree: 


i=l 


+25. (p 2)) er) yP 6) (a - >on) (2.20) 


where 6,(p?) = 6(p?)@(po) with @(po) being the Heaviside step function which imposes 
the positivity condition on po. In appendix E, we give explicit parametrizations of the 
two-particle and three-particle phase-space integrals. 


“Note that in (2.19) the complex conjugate of tree-level form factors is already encoded in replacing O 
by its conjugate O and changing the MHV degree from k to m — k. Therefore, the ‘*’ refers to taking the 
conjugate of the £ > 1 contributions only. This will be explained in explicit examples in section 5; see the 
discussion around (5.6). 
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Finally, the sum over the different channels, i.e. over the different particle numbers n, 
has to be performed. This leads to a cancellation among the different soft and collinear IR 
divergences such that the final result is IR finite. If non-protected operators are involved, 
as in the Konishi case, their renormalization constants have to be taken into account. 


3 Form factors for Ke via unitarity 


In the previous section, we have defined the cross section for gauge-invariant operators O 
in MN = 4 SYM theory in terms of its squared matrix elements. As discussed around (2.19), 
the building blocks of these squared matrix elements are the non-color-ordered super form 
factors for the respective operator. In this section, we will present the building blocks 
necessary for computing the cross section of the Konishi operator (1.4) up to two loops, 
which are the two-point form factor up to two-loop order and the three-point form factor 
at one-loop order.!? 

We use the notation FS) (Lai, pi; m}) for the non-color-ordered super form factors 


and FY ({pi, mid) for the color-ordered super form factors, as introduced in (2.11). We 


denote the bosonic color-ordered form factors with fixed external states by Fe ) ({pi}). If 
necessary, we specify the external states by subscripts, e.g. in case of two scalars and one 


gluon as Fe Hig 24,39) or simply F a These bosonic form factors can be obtained 


b,0,9)° 
from FY Upi mY) by taking a specific term in the 7; expansion. We also introduce the 
normalized bosonic form factors E ) as the ratio between the 4-loop and tree-level color- 
ordered bosonic form factors: 


(£) 
(0) Ay Fo (pis) 
fòn Pi) FO (p) $ (3.1) 


Our computation will focus on the colored-ordered form factors; via (2.11), it is straight- 
forward to obtain the full non-color-ordered super form factor from them. 

The computation of form factors in this section are based on the on-shell superspace 
formulation (1.7). Therefore, the operator in the form factor is Ke defined in (1.10) and 


not the Konishi operator K defined in (1.4). We denote the resulting form factors by F oe 


) 


7» we have to modify the results presented 


In order to obtain the Konishi form factors Fe 
in this section, as will be discussed in detail in section 4. 


3.1 Some BPS form factor results 


We start by presenting some known results for BPS form factors, which are also useful 
building blocks for the Konishi form factors. Unless otherwise specified, the BPS form 
factor in this paper will always refer to that of the half-BPS operator tr(¢4, pg) defined in 


(1.9), and we use the abbreviation FE, — (8) 5 : 
’ tr(¢4_).n 


"The tree-level four-point Konishi form factor essentially agrees with the BPS result, as we will discuss 
in subsection 3.2. 
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The n-point MHV tree-level BPS super form factor is given by [22] 


sO aa) 
W223). ml)’ 


(0), MHV 


Fpps.n (1,2; a nN) = —_ (3.2) 


where 6(4)4B (X; Aini) is understood as taking 77 in the delta function with only A, B indices, 
or more explicitly 


wae S dum) = (SG ann) (eon?) > (3.3) 


i<j k<l 


Note that in this and all following expressions for form factors we do not explicitly write the 
momentum-conserving delta function (27)*6“) (q— X}; pi), where q is the four-momentum 
carried by the gauge-invariant operator. 

We give the loop corrections to the BPS MHV form factor in terms of the normalized 
form factor defined in (3.1). In this paper, we only need the following three results [13, 20, 
21): 


= —2819 (3.4) 


i 2- 


E 2 = Sip (3.5) 


pı 

1 ee _ $13 + $23 5 523 _ $12 $ $12 + 531 
a=- E Taci <Ç 
P3 


+ cyclic perm. of {p1, p2, P3} , (3.6) 


where sij...k = (pi + pj +: + pk}. Each graph corresponds to a Feynman integral which 
is defined in appendix B. Throughout this paper, all external on-shell momenta p; are 
understood as outgoing. 

For the two-point case only the MHV configuration exists, while at three points there 
are the MHV and the next-to-MHV (NMHV) configuration. The NMHV tree-level form 
factor can be obtained from (3.2) by first taking the conjugation \ > \ and n4 > ña, and 


then applying a fermionic Fourier transformation ast? 


- ne 4b jet Mit) 
(0), NMHV _ de Ch: \ TAB\Cg=1 ^ii 
Fon 1,2,3)={ ]| pek tio |S ; 
me (2,3) (Ife o ) n231] an 


The loop correction to both, the MHV and the NHMV three-point form factor, is given by 
(3.6). 


13Recall that the operator also becomes the conjugate one, tr((¢47)?), where 47 = ze4BCP bop. 


Ee 


3.2 Tree-level two- and three-point form factors 


We now turn to the form factors of Kg. In this subsection, we consider its tree-level form 
factors. They are identical to those of the Konishi operator K. The expression for Ke 
in (1.10) contains the individual fields ¢4g¢cp where A, B,C,D assume distinct values 
instead of ¢łp as is the case for the BPS operator. For the tree-level bosonic form factor 
with specified external particles, however, the index structure of the external scalars and 
fermions do not play any role in the result, which is obvious from the Feynman diagram 
computation. Therefore, the tree-level bosonic form factors for the Konishi operator are 
identical to the corresponding BPS form factors. 

The super form factors, on the other hand, take different forms. Taking into account 
all the components, the two-point super form factor reads!* 
1 12} 


FR (2) = -jaganing ne?) > (3.8) 


where £1234 = 1. The bosonic two-point form factor 


2 
Be Dia) = ar =1 (3.9) 


can be obtained by taking the (n{n?)(n3n$) component of the tree-level form factor Fs 
in (3.8); it is identical to the BPS result as can be seen by taking the (nin?) (én?) 
component of (3.2) at n = 2. There are two other possible scalar field configurations at the 
external legs, namely {(¢13, 24), (14, ¢23)}, and for both these cases we obtain the same 
bosonic form factor as above. 

The three-point MHV super form factor is given by the following expression: 


FE) 2.3) = grygy (l 2 ason nin!) (ine) a 


+ 2(13)(23)eancpntink (n ny) + cyclic perm.) 


It has two distinct configurations of the external states: scalar-scalar-gluon and fermion- 
fermion-scalar. Taking the coefficients of (n{n?)(n3n3) and nin (nèng), we find 


FO (14,5, 24043 a o POG. E a 
Ke \~o127 A347 99+ (12)(23)(31) ’ Ke \ Y1? 29 Sosa (12)(23)(31) ’ 


(3.11) 
which are also identical to the corresponding BPS form factors. The NMHV form factor 
can be obtained from the MHV result in a similar way as in the BPS case (3.7). 


3.3 One-loop two-point form factor 


In this and the following subsection, we compute the form factor of Kg at one- and two-loop 
level via four-dimensional unitarity [3, 4]. 


“The normalization factor is fixed to be consistent with the definition of the operator Ke in (1.10). 
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Pı 


P2 
Figure 2: The simple (pı + p2)? double cut. 


The general idea of unitarity in this context is to reconstruct loop corrections to the 
form factors at the integrand level from their discontinuities, i.e. by applying cuts. Here, 
a cut denotes setting a propagator on-shell according to 


2 + 2784 (l2) , (3.12) 
where 6,(1?) was defined after (2.20). On the cut, the loop expression factorizes into a 
product of (known) tree-level or lower-loop form factors and amplitudes. These have to be 
summed over all possible particles exchanged in the cut channel, which can be achieved by 
integrating the super form factors as well as the super amplitudes over the Grassmannian 
degrees of freedom in the cut legs. Then, one can apply the spinor algebra to write the 
result in a form that can be identified as a sum of cut integrals. In this way, an ansatz 
for the uncut integrals occurring in the loop correction is assembled. In general, not all 
integrals appear in a given cut, and additional cuts have to be taken to complement the 
ansatz. The complete ansatz has to be consistent with all possible cut. Finally, the cut 
integrals have to be lifted to the uncut integrals, as discussed in appendix B. 

In the following, we apply this technique to the form factor of Kg. We start with the 
computation of the one-loop two-point form factor. For the sake of explicitness, we choose 
a fixed combination of external scalar states, namely {¢12, 434}. As in the tree-level case, 
the other two choices of external scalars {¢13, 624} and {¢14, 623} lead to the same result. 


We abbreviate es 2434) as ph (6oy 


Only one cut needs to be considered: the two-particle cut in the channel (pı + p2)? = 
q?.15 It cuts the internal propagators carrying momenta lı and lọ as shown in figure 2. 
The building blocks on the two sides of the cut are the color-ordered two-point form factor 
(3.8) and the color-ordered four-point MHV amplitude given in the standard MHV form 
[68] as*® 


SOO Am) 
An =T aI an) 


The other two two-particle cuts occur in the p? and p3 channels. Since these legs have p? = p3 = 0, 
massless bubble integrals in these channels vanish identically when regularized in D = 4 — 2e dimensions. 
Hence, all integrals can be detected by the q?-cut. 

16Recall that we are always suppressing the momentum-conserving delta function in the notation. 
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The sum over all possible particles exchanged along the cut is considered by integrating 
over the fermionic coordinates of the exchanged particles as f dôn, , = f dm, dm, while 
keeping the external state fixed. 

The q?-cut integral reads!” 


= | dPSo47 Bm, Fo (h, -l2) x AP (p1, pa, lo, 
: f 2,1} 4 Maı,2 Ko 2k 1, —l2) 4 (Pi: P2, l2 Deg 17 )(n32) 


(l1 2)? (l2 1)? + 4b 1) (l1 2) (l2 1) (l2 2) + (l1 1)? (l2 2)? 
=V , yi / dPSo, pH FEU(CEIC TS e 


C(4,0) 


(1) 
PKG, ’ 


(3.14) 
Since the external states are fixed to be {d12, 631}, we take the (n}n?)(n3n3) component 
of the cut integrand. The phase-space integration measure, dPS2 7}, is defined according 
o (2.20), with the integration variables being the momenta of the cut propagators {l1, l2}; 
hence the subscript in the notation for dPS 7. 
The cut integral can be simplified at the integrand level as!® 


_ (lila)(12) (lv 2) (la I) \ _ — s12 (l + po)? 
Cog) = J aPS2.(0 G (a2) Aa) 5) ~ / i (z FoP sip ) 
l i 


1 pı 5 i pı 
D << paltar X 
i. p2 P12 k p2 

Š (3.15) 
where the flow of the momenta is as specified in figure 2. In the above equation, the 
integral over the two-particle phase space is shown by the dashed cut line of the triangle 
and bubble graph. For the triangle graph, the denominator in the integrand is the uncut 
propagator Gin? and the numerator coefficient is —sy2. The shown bubble graph has 
no uncut propagator, but is has a loop-momentum-dependent numerator factor, which is 

written in front of the graph. 
As described in appendix B, the cut integrals (3.15) can be lifted to the full integrals. 


The full normalized form factor as defined in (3.1) then becomes!’ 


pı l pı 
Do o = sa 
keeg = 2 | 752 , ra X, ' (PA 
2 2 


where the factor of 2 is due to the permutation of the two external legs, and we use the 
short notation sj = (pi + 1). Note that the prefactors that depend on the loop momentum 
are understood to appear in the integrand of the integral represented by the respective 
graph it multiplies. 


'7For reversed momenta | — —l, occurring e.g. in i ) 2(- lı, —l2), we follow the convention A; > —A; 
and Ài > Ni in the spinor helicity formalism. 

®The first line can be obtained via the Schouten identity (ab) (cd) = (ac) (bd) + (ad) (cb) for (l 2) (l2 1) 
n (3.14). 

19The coupling dependence can be recovered as shown in appendix B. 
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The bubble integral with loop momentum in the numerator can be reduced to the 
scalar bubble integral via Passarino-Veltman (PV) reduction, see appendix C for details. 
Thus, we obtain the form factor? 


Pı PL 
fe = Stay <i 6 -OX l 
Ke,( $,0) 7 bs (3.17) 


=. — 
ee 


The integrals corresponding to the graphs are given in appendix B. Note that the contri- 
bution to the form factor involving the triangle integral is the same as the BPS form factor 
23 in (3.4). An independent computation of this result via Feynman diagrams is shown 
in appendix H. 

From the above calculation at one-loop, we see that the IR-divergent part of the form 
factor of Kg is the same as the one of the BPS operator. The extra contribution coming 
from the UV divergent bubble integral yields a non-vanishing anomalous dimension unlike 
in the BPS case. We will equally organize all subsequent results for the form factor in 
terms of a part that is identical to the BPS form factor and an additional contribution that 
is unique to the form factor of Ke. 


Vanishing one-loop form factors 


Before proceeding to two-loop order, we briefly discuss two other possible form factors with 
gluon or fermion external states, namely the form factors F D (lg_,2g,) and F, 3 (Lyri 2z): 
At tree level, they are zero since no Feynman diagram for this configuration exists. At 
higher loops this is not obvious. Here, we use unitarity to show explicitly that they are 
zero at least at one-loop order. Consider the q?-cut as in (3.14), but to obtain ia (la; Zax) 
and TARER take the components nin?n?nf and nt(n3n3nz) of the cut integrand, 
respectively. This yields for the two cases 


ea (Aa? aor 
, =i f aPs, ET FE <{ , (3.18) 


„=i f aps 3(l1 1) (l1 2) (l2 1)? + 3(l1 1)? (l2 1) (l2 2) 
ae (12)(h o) 1) (2) 


lı pı i pı 
' 1412 

=i3(1|4|2] =< pag loll -0X l (3.19) 
i p2 $12 ‘ p2 


2 le 


When we lift these expressions to the full triangle and bubble integrals and perform the 


Feet 9234) 


PV reduction, we obtain zero. Since we use four-dimensional unitarity, we also have to 
check that there is no contribution from potential rational terms. A similar (but simpler) 
study as in appendix D shows that rational terms are indeed absent. 

Finally, there is an easy way to see that Ficg(1g_,2,,) = 0 to all loop orders. Using 
the gauge freedom, we can choose the polarization vectors of the outgoing gluons as £] = 


?0For convenience, we will from now on refer to the normalized form factor as form factor, too. 
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eg x À dg. It is then obvious that the form factor must be zero, since it is proportional to 
Ei * Pj OF E1 * E2. 


One can also compute FY 


Kela) directly by using Feynman diagrams. A simple com- 


putation gives 


E1 - p2) (E2 p 
Fup = 2e ea) ERE pe, (3.20) 


where the integral I? [Z] = 5 + O(e) is given in (B.10) for pp —> 0 and the relabeling 
p3 > p2. This result holds for the polarization vectors ey taken to be in general D = 4— 2e 


dimensions. Since J3[¢?] is finite and its prefactor is of order O(c) (as it vanishes when 


D = 4), the form factor itself FO 


«.(9.g) 18 of order O(e). This is consistent with the unitarity- 


based calculation. 


3.4 Two-loop two-point form factor 


Next, we compute the two-loop two-point form factor of Kg. As in the one-loop case, we 
specify the external states to be {¢12, d34}. 


Two-particle cut 


We first study the two-particle cut in the q?-channel. We follow a similar procedure as the 
one being used in computing the BPS form factor [27]. We first quote the q?-cut integral 
given in (2.6) of [27]: 


FSM) = J PS2 dm FS h fa) (4AP (Pr pas last) + AY Prl p2:b)) » 
(3.21) 
where the building blocks are the two-point tree-level form factor (3.8) and the one-loop 
color-ordered four-point amplitude [69] 


AP (p1, p2, p3, pa) = AY) (p1, P2, P3, Pa)(— 512823) I$} (p1, po, p3, pa) : (3.22) 


The tree-level super amplitude AW) (p1, P2, p3, pa) in (3.22) contains all the dependence on 
the fermionic coordinates, and the term multiplying it is a massless scalar box integral [ 3 
defined in (B.9).?? 


Let us briefly explain (3.21); see [27] for a derivation in full details. The above cut 
integral is obtained by taking the product of the two-point form factor and the non-color- 
ordered four-point amplitude. The one-loop four-point amplitude contains a single-trace 
contribution, as well as a double-trace contribution which is sub-leading in color. However, 
after the contraction of the color factors with the two-point form factor, both contribute 
to the cut integral with the single-trace color factor ô% = tr(T*T°).”° The final building 
blocks in the cut integral are the color-ordered form factor and amplitude as given in 
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Pı 


P2 
Figure 3: The (pı + p2)? double cut at two loops that contributes to the planar ladder 


integral. The building blocks are the color-ordered tree-level form factor and the color- 
ordered one-loop amplitude. 


(3.21). The two contributions in the parentheses of (3.21) are depicted in figures 3 and 4 
respectively.2+ We consider them one by one below. 
We first study the contribution from the first term in the parentheses of (3.21), which 
La is shown in figure 3. Using the one-loop result (3.22) for the amplitude and taking the 
> external states to be {¢ġ12, $34}, the corresponding cut integral can be written as 
Do 8 (0) (0) 
FE 6.) | T f Psa d mmaFKe2(lh la) Ag (pi Pasta, 11) (nin?) (nn) (3.23) 
x (—)s12811, 1)? (p1, p2; l2,l1) - 


The first line in (3.23) given by the product of tree factors is the same as the cut integrand 
of the previously studied one-loop case in (3.14). Therefore, we can perform exactly the 
same calculation as in (3.15) and obtain 

I 


(2) _ (0) —$12 $21 (1) 
FG (,0) e = =F kele) f Psa S11, + 6) sias, 14 (p1, p2, l2, l1) 


li pı 
= FK (4) (Sta — 6s san ) | 
l p2 


2 


(3.24) 


The above cut integral can be lifted to the two-loop planar ladder integral. This 


integral can be drawn in two different ways, namely 


pı pı 
<1 , Sl] (3.25) 
p2 p2 


21Note that (3.21) applies to any composite operator with two elementary fields, in particular to Ke. 

22 The minus sign in (3.22) is related to the convention of the box integral we use in (B.9). 

?3The enhancement of the power in Ne of the apparently suppressed double-trace term in the amplitude 
is the wrapping effect analyzed earlier for the spectral problem [70]. 


24The factor 4 in the first term comes from the different contributions of the color factor contraction; two 
of them come from the single-trace four-point amplitudes, the other two from the double-trace four-point 


amplitude, as explained in [27]. A different way to understand the factor 4 is to look at the two-particle 
cut with the one-loop form factor on the left hand side and the tree-level amplitude on the right hand side. 
It then arises from twice applying the reasoning that gave us the factor 2 at one loop. 
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Figure 4: The q?-cut at two loops that contributes to the crossed ladder integral. The 
building blocks are the color-ordered tree-level form factor and the color-ordered one-loop 
amplitude. 


Furthermore, there are two other planar graphs obtained by permuting the external legs 
pı © p2. So, altogether we have four diagrams which are drawn in different ways but 
all give equivalent planar ladder integrals. This provides a diagrammatic interpretation 
of the factor 4 in the first term of (3.21). As we will see later in the triple cut, it is also 
important to separately draw the ladder graphs in different ways according to (3.25) in 
order to compute the cut integrand correctly. 

Next, we consider the second term inside the parentheses in (3.21), which is depicted 
in figure 4. The corresponding cut integral is given by 


(2) 
PK old, al. 


=- f aPs,, ty dên, Popa l, ~1) AO (pı, 11, pa, l2) 
(3.26) 


(1) 
X 8198 I „l ’ l : 
128114 (Pı, 41, pa, la) (nin? )(n3n3) 


Following similar steps as in the previous case, the cut integral is expressed as a two-particle 
cut of the two-loop crossed ladder integral as shown below 


2 II 0 
’ Ke.) ep = FS 66) (si. = 6s san ) P (3.27) 


Now, lifting the cut integrals of the combined contributions (3.27) and (3.24) as described 
in appendix B, we find the following contribution to the two-loop form factor of Ke: 


l Pı l pı 
(2),I (2), II 
Eo + Sega) = (82 — 6susa) eT] + << , (828) 
p2 p2 


where, as explained around (3.25), the factor 4 is included for i dy" 
There is another g?-cut which is similar to the one in figure 3. It has the one-loop 
two-point form factor on the left hand side and the tree-level four-point amplitude on the 
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Pı 


P2 
Figure 5: The two loop (pı + p2)? = q? triple cut. 


right hand side. This case is a bit subtle. Naively, one would expect that only the one- 
loop two-point form factor with scalar external states Fẹ (1) Ke,(6,0) CaN occur on the left hand 
side, since the other possibly contributing form factors with gluon and fermion external 
states FU e PA 
However, this expectation turns out to be incorrect. There are non-vanishing contributions 


(1) ; f 
and F Ke (h1,p234)? , respectively, vanish as shown in the previous subsection. 


from these two cases: only the integrated one-loop form factors are zero, but the integrands 
are not, as we can see from (3.18) and (3.19). Their integrands have to be taken into account 
in the unitarity cuts and then yield a result which is consistent with the one found from 
the q?-cut of figure 3.” 

The result (3.28) obtained by using only the two-particle cuts is not guaranteed to give 
the full form factor. One problem is that the numerator coefficients of both integrals are 
ambiguous w.r.t. terms that are proportional to /?. Due to the on-shell condition of the 
cut propagators, such terms are not detected by the double cuts. Moreover, there may be 
other basis integrals which cannot be detected by the double cuts. Both these issues can 
be fixed by studying the three-particle cuts, which we do next. 


Three-particle cut 


The three-particle cut, or triple cut (TC), across the q?-channel is shown in figure 5. Unlike 
for the BPS form factor, the triple cut will indeed give some new contribution to the form 
factor of Kg, which is not detectable by the previous double cut. 

The cut integral is given as 


), MHV .NMHV 
FO otha oo [PSs ni mF WY (h, =z, =l) AP (p1, p2, l3, l2,l1) 
Fom In, =la, =s) AP EY (p1, pa, ts,120)) | ane 
1/1 2°12 
(0) 
te ($:$) WUTA (3.29) 


Note that besides the MHV form factors and amplitudes also the NMHV form factors and 
amplitudes appear as building blocks in (3.29). The two terms in the above sum are in 
fact conjugate to each other. 


?5Since the non-planar ladder does not contribute to this cut, one only obtains the contribution coming 
from the planar ladder integral in (3.28). 
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Figure 6: Triple cut of the integrals that correspond to the first five terms in (3.30). The 
flow of the momenta is as specified in figure 5. 


After performing the fermionic integrations? and some spinor algebra, the cut integral 


can be simplified at the integrand level to obtain the following form:?’ 


2 2 2 
Cro = [ass m (= = 6811, $21, + $19 — 68113 5213 + S12 — 68115 S212 
; $213 Sli lz Slal3 $11, Sly lz Sləl3 $213 Slı lə Sli l3 


87, — 6811821, n 845 — 6811p S212 x 18 18s, — 18s2, ) 
$11, Slılz Sləl3 $11, $213 Slıl3 $12 8128lylg  $12Sləl3 
(3.30) 
Note that the first five terms in (3.30) can be obtained directly from the result de- 
termined by the two-particle cut in the previous paragraph, namely the planar ladder 
contribution in (3.24) and the crossed ladder in (3.26). Let us first look at the first term in 
(3.28), the contribution from the planar ladder, which contains the numerical prefactor 4. 
As mentioned earlier, this factor 4 stems from the four different ways of drawing the planar 
ladder graph. The two configurations shown in (3.25) contribute to the above triple cut. 
In order to account for all possible triple cuts on these two diagrams, we cut each in two 
ways as shown in figure 6. Thus, the first four terms in (3.30) correspond to the first four 
diagrams in figure 6, which are just the planar ladder integrals. The remaining fifth term 
in (3.30) correspond to the last diagram in figure 6, which is the crossed ladder integral 
with only one possible triple cut. Hence, the first five terms in (3.30) do not result in any 
new contribution but reproduce the double-cut result in (3.28). 


The remaining three terms in (3.30), however, are new contributions to the two-loop 
ansatz detected by the three-particle cut. They can be expressed as the three-particle cut 


26To obtain the cut integrand in a compact form, it is convenient to take the product of the bosonic 
form factor and amplitude expressions and sum over all helicity configurations, since the NMHV result of 
both, the three-point form factors and the five-point amplitudes, take simple MHV form. We have checked 
that the expression obtained in this way is equivalent to the expression by using super form factor and 
amplitudes and doing the fermionic integration directly. 

27Tn practice, this form can be obtained easily as follows. First, one can write down immediately the 
contribution of the first five terms by using the result (3.28) obtained from the double cuts, as explained 
below. Then, subtracting them from the cut integrand, the remaining terms take a very simple form which 
can be easily simplified into the last three terms. 
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These three cut integrals can be lifted to full integrals, which can be simplified further at 
the integral level to give a single scalar integral: 


pı 
2), II 
a -87X ] (3.32) 
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Complete two-loop result 


Now, we combine the results from all the cuts, (3.28) and (3.32), and obtain the two-loop 
two-point form factor,” 


E (2),11 (2), 11 
ape ee are ae 
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(3.33) 


where the integrals corresponding to the graphs are given in appendix B. Note that we have 
multiplied f a 7 
legs pı © po. wee in the one-loop case, we have presented the result by separating a part 
that is identical to the BPS form factor fos given in (3.5). 

The double and triple cuts we have considered should be able to detect all possible basis 


$) by 2 to include the contribution from the permutation of the external 


integrals up to potential rational terms that might be missing when using four-dimensional 
unitarity. Comparing our result (3.33) with the one we obtained for Ke from the Feynman 
diagrams of appendix H, we have confirmed that such rational terms are absent. 

As will be explained in section 4, the result given by (3.33) is, however, only valid for 
the operator Ke defined in (1.10), but not for the Konishi operator K defined in (1.4). This 
subtlety will be discussed in details in section 4. We will see that by a rigorous prescription 
we can modify the above result in order to obtain the Konishi form factor. 


3.5 One-loop three-point form factor 


In this subsection, we compute the one-loop three-point form factor of Kg. The computation 
is similar to what we have done for the previous two-point case. We need to consider cuts 
in all possible kinematic channels, which, apart from the q?-cuts employed earlier for the 


?8 The result (3.33) matches the one in the unpublished notes of Boucher-Veronneau, Dixon and Penning- 
ton [71]. 
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Figure 7: The cuts needed to compute the one-loop three-point form factor of Ke. 


two-point form factors, contain also the sab-cuts, where (a, b) = (1, 2), (2,3), (3, 1), as shown 
in figure 7. Combining the results from both types of cuts ensures that no contribution to 
the ansatz is missed. 


Unlike for the BPS form factor, the loop corrections of the form factors of Kg turn 
out to be different for different configurations of external particles. Therefore, we need to 
consider the form factors with specific configurations of the external states individually. 
We consider the scalar-scalar-gluon and fermion-fermion-scalar cases.7? We will discuss the 
scalar-scalar-gluon case in some detail. The fermion-fermion-scalar result can be obtained 
in the same way and we only present the final result. 


1 
FE) gaz , 2¢cp 9 394 ) 


We first consider the form factor of Kg with scalar-scalar-gluon external states. For the 
sake of explicitness, we focus on Feces Oe ,), which we abbreviate as oe 4-b.0)" 
The result applies to all other non-vanishing cases, where A, B,C, D are distinct and the 
g may have positive or negative helicity. As shown in figure 7, we need to consider both 
the q?-cut and the s,y-cut. Since the operator is a color singlet, we need to consider all 
possible cyclic permutations of external on-shell legs in the cuts, as they contribute to the 
same color-ordered form factor. Explicitly, we need to consider three cases for each channel 


in figure 7: 


{a, b, c} = (1) 1 Libis Abana 394) (II) {2345 394» léit, (IIT) {39,, Liss Zost f (3.34) 


In total, there are six cut channels to consider: (a-I), (a-II), (a-III) and (b-I), (b-ID), (b-III), 
where (a-I)-(a-III) are the q?-cuts while (b-I)-(b-IIT) are the s,)-cuts. Note the (I) and 
(III) cases are actually related to each other by a flipping symmetry. 


?°There could also be other external states composed of different fields, such as F de ,29,,39,) and 


PO) iis Bii 3g, ). Such form factors, however, do not contribute to the two-loop cross section studied 
in section 5 as the corresponding tree-level results are zero, and we will not consider them in this paper. 
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(a-I)-cut: 


This is the q?-cut in figure 7 with the choice of external legs {pq, pp, Pe} corresponding to 
the particles {14,,,243,,39,}- The cut integral is given by the following equation: 


(1) 
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C(a-1) 
(3.35) 
where the tree-level form factor res CC rps ee 2434139, ) is given in (3.11). 
The above result can be reduced to an appropriate cut of integrals by using some 


spinor algebra. Without going through the detail, we present the result: 
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This is the q?-cut in figure 7 with the choice of external legs {pq, Pb, pe} corresponding to a 
different order of particles, namely {24,,,3,,,1¢,,}. The cut integral can be computed as 


(a-II)-cut: 
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(3.37) 
After some spinor algebra, the above result can be expressed as cut of the following inte- 
grals: 
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(b-I) cut: 


This is the sı2-cut in figure 7 with the choice of external legs {pa, pp, Pe} corresponding 
to the particles {14,,,2¢43,,39,}. In this case, one of the building blocks is the tree-level 
three-point form factor in (3.10). The cut integral is given by, 


1 0),MHV 0 
TE E = f PS den, Fe (— hh, Ip, 3) AS ) (p1, po, lo, l) 


(nin?) (nn) 
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(3.39) 


It can be written as cut of integrals as 
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This is the case of the sg3-cut, the last of the independent cut channels, in figure 7 
with the choice of external legs {pa,py,p-} corresponding to the particles in the order 
of {2434,39,,1¢,.}. The cut integral is given by 


(b-II) cut: 
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As previously mentioned, the cuts (a-III) and (b-III) give similar results to (a-I) and 
(b-I) and can be obtained by exchanging pı + pə in the latter. Hence, we will not give 
them explicitly. 
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From the cuts to the full form factor 


We find that all the above cut results are consistent with each other, i.e. the prefactors 
of the graphs are identical when the same graph appears in different cut channels apart 
from terms that vanish due to the on-shell condition for the cut momenta in the individual 
channels. Given all these cut results, it is straightforward to lift the cut integrals, as 
described in appendix B, to obtain the full form factor®? 


pı 
a) 7 (s123 — 2513) (s13 + $11) 13 (813831 — S128521) 1 
f = ; 3| > - p2 
K6,(¢,9,9) a 5235123 2 p3 
3891 (s11(513 + 523) + jo : 
$91 (S11 (813 + $23) + 512523 l 
ee eee h +p o pa} Aa 
12 p3 


(3.43) 
where f a denotes the BPS part that is given in (3.5). 
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For the form factor with fermion-fermion-scalar external states with distinct A, B, C, D, e.g. 
FY 
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6 
in all possible channels and lifting the cut result to the full answer. Without giving details, 


(Lyi 242) 3034), One can proceed along the above steps for computing the cut integrand 


we present the final result, denoted by oa abg) 
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Note that in the above result not all contributions from box graphs are incorporated in the 


(3.44) 


corresponding BPS part given in (3.5), unlike in the expression for the scalar-scalar-gluon 
form factor (3.43). Moreover, there is an additional one-mass triangle integral in the first 
line of (3.44), which does not appear in (3.43). 


PV reduction and some interesting features of the results 


We have obtained the full integral expressions for the form factors (3.43) and (3.44) of Ke. 
The results are obtained by using the unitarity method fully at the integrand level. As 


30 Recall that the prefactors that depend on the loop momenta are understood to appear in the integrand 
of the integral represented by the respective graph each prefactor multiplies. 
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a result, the integrals still contain loop-momentum-dependent numerators. Such integrals 
can be reduced further via PV reduction, see appendix C for details. 
After PV reduction, the results (3.43) and (3.44) are simplified to 
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where all relevant integrals are given in appendix B, and Fin extracts the finite part FB of 
the box integral defined in (4.8). 
There are several interesting features in the above results we would like to comment 


AES Soe (3.45) 


on. 


e Going back to the expressions (3.43) and (3.44), we notice that — besides a part 
identical to the BPS form factor a — there are still triangle or box integrals left, 
which separately contain IR divergences. This might cause a net IR divergences in 
addition to the one contained in o i.e. it would spoil the universality of the IR 
divergence. However, as evident from the results (3.45) and (3.46), these additional 
IR divergences cancel after PV reduction and hence only the universal IR, divergence 
of the BPS part remains; see also [20]. 


e There are remaining divergences given by bubble integrals. These are the UV diver- 
gences which have to be canceled by renormalizing the composite operator Kg. See 
section 5 for a further discussion. 


e Besides the common BPS part, the two results (3.45) and (3.46) are quite different 
from each other. This directly shows that the form factors of Kg with different 
external legs (even with the same MHV degree) have very different structure and 
need to be studied case by case. 


e There is a term in (3.43) involving the integral IP [¢2] given in (B.10). It evaluates 
to a rational term. Interestingly, we have found this contribution by applying four- 
dimensional unitarity. This is possible, since we apply the four-dimensional unitarity 
to compute the integrand expressed in terms of a tensor-integral basis. The rational 
term only appears after the PV reduction when the basis is reduced to scalar integrals. 
In the usual one-loop (generalized) unitarity computation [72, 73], one computes the 
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coefficients of the scalar integrals directly and hence one would miss this rational 
term. In appendix D, we have checked that the rational term in (3.43) matches with 
an independent Feynman diagrammatic computation, and hence the final result for 
the form factor is complete. 


e The integral coefficients of the form factors in (3.45) and (3.46) appear to contain 
EEE = ae These are, however, just spurious poles 


which cancel when all contributions are taken into account. We demonstrate this in 


unphysical poles, such as 


detail in appendix D. 


Finally, we would like to mention that most of the above features (except the IR divergence) 
do not occur for the one-loop scattering amplitudes and BPS form factors of the NV = 4 
SYM theory. In QCD, they are, however, common and appear e.g. for one-loop amplitudes 
[5]. 

Last but not least, recall that the above form factor results for Ke still need to be 
modified to obtain the correct ones for the Konishi operator, as will be described in the 
next section. This does not affect any of the above listed properties. 


4 Konishi vs. Ke 


In this section, we discuss some important subtleties that arise when regulating the theory 
by continuing the spacetime dimension from D = 4 to D = 4 — 2e. Our unitarity-based 
calculation made use of the on-shell superfield (1.7) that only captures all degrees of freedom 
in strictly D = 4 dimensions. Hence, this approach does not directly yield the correct form 
factors of dimension-dependent operators. We explain this problem and its resolution in 
details below, taking the Konishi form factor as a concrete example. 


4.1 A subtlety in choosing a regularization scheme 


When regulating the theory by continuing the spacetime dimension to D = 4—2e, one has to 
also specify how the various fields are continued. In conventional dimensional regularization 
(CDR) [74] and the ’t Hooft Veltman (HV) scheme [75], the number of fermion flavors Ny 
and also the number of scalar flavors Ng remain as in four dimensions and are hence kept 
as Ny = 4 and Ng = 6, respectively. This does, however, break supersymmetry, since the 
polarization vector é“ is taken in D = 4 — 2e dimensions. 

A scheme that preserves supersymmetry is dimensional reduction (DR) from ten di- 
mensions [42, 43]. In this scheme, the number of scalar fields is changed to Ng = 6+2e, such 
that D + Ng = 10 is independent of e. It exploits the fact that four-dimensional M = 4 
SYM theory can be obtained by dimensional reduction of ten-dimensional M = 1 SYM 
theory. Performing the dimensional reduction to D = 4 — 2e rather than four dimensions, 
one obtains a regulated theory that preserves V = 4 supersymmetry. The ten-dimensional 
gauge field Am, M = 1,...,10, then reduces to the D-dimensional gauge field A, and to 


MN 


Ng = 10 — D = 6 + 2e scalar fields øz. Similarly, the ten-dimensional metric g reduces 


to the D-dimensional metric g” and 6/7. 
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In a modified version of the DR scheme, known as four-dimensional-helicity (FDH) 
scheme [76, 77], the 2e scalar degrees of freedom are absorbed into the gluons. This 
apparently preserves supersymmetry in the sense that bosonic and fermionic degrees of 
freedom still match. In particular, it allows to use the M = 4 on-shell superfield (1.7) and 
polarization vectors in D = 4 dimensions. So far, the FDH scheme has been successfully 
used in computing amplitudes and BPS form factors in M = 4 SYM theory. However, as we 
will discuss below, the FDH scheme is incompatible with dimension-dependent operators, 
i.e. operators that are sensitive to the absorption of the 2e scalar degrees of freedom into the 
gluons. In particular, the Konishi operator (1.4) is dimension-dependent, since it contains 
a trace over Ny = 6 + 2e scalars in D = 4 — 2e dimensions. The incompatibility arises 
in the FDH scheme since the 2e scalars are absorbed into the gluons, giving direct results 
only for Ke in (1.10). 

In order to detect the differences between working with the FDH and the DR scheme, 
we examine the underlying Feynman diagrams. In Feynman diagrams, explicit factors of 
D = g’, and Ng = ôl arise whenever a gauge or scalar field runs in a loop in which the 
respective Lorentz or flavor index also forms a loop. We call such a loop an index loop. 
Moreover, we call an index loop externally closed if fields of the composite operator are 
involved in the index loop and internally closed if the operator is not involved in the index 
loop.?! 

We look at internally closed index loops first. From the dimensional reduction from ten 
dimensions, we know that an internally closed vector index loop always occurs together 
with an internally closed scalar index loop. This is illustrated in figure 8. Hence, each 


Figure 8: In dimensional reduction, a Feynman diagram with a closed ten-dimensional 
vector loop (zigzag) decomposes to one with a closed 4(—2e)-dimensional vector loop (wig- 
gly) and one with a 6(+2e)-dimensional scalar loop (plain). 


factor of D for a Lorentz index loop is accompanied by a factor of Ng for a scalar flavor 
index loop. This can also be seen in the concrete Feynman diagrams in appendix H, e.g. by 
comparing the first two lines in table 2. The sum of both contributions is proportional to 
D + Ng = 10, both in the DR scheme and in the FDH scheme. Hence, as far as internally 
closed index loops are concerned, one is free to work in the FDH scheme. 

The situation changes for externally closed index loops. Generically, the fields of a 
composite operator involved in such a loop are only a subset of the fields in the theory, e.g. 
only the scalar fields. In this case, a diagram in which the externally closed scalar loop 
generates a factor Ny is not paired with a diagram in which a vector field can circulate in 

31Note that the factors of D = g”, and Ng = 67; occur even if the index loop is apparently interrupted 
when the gauge or scalar field splits into a pair of fermions that themselves build a loop. This follows 


from the Clifford algebra for the spacetime (o,)aq matrices and the flavor (07) 4B matrices into which the 
ten-dimensional ø matrices split. 
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the loop and generate a factor D. Hence, the result in the FDH scheme differs from the 
one in the DR scheme. In scattering amplitudes and form factors of BPS operators such 
as tr(¢*,), no externally closed index loops can occur. This is why the FDH scheme is 
directly applicable for calculating these quantities. 

Let us consider the particular example of the Konishi primary operator (1.4); it is 
defined as the trace over all scalars and can hence be part of an externally closed scalar index 
loop. The Konishi primary is the highest-weight state of a so-called long supermultiplet of 
psu(2, 2/4). Supersymmetry guarantees that all members of this supermultiplet have the 
same anomalous dimension (1.6) — unless it is broken by the regularization scheme. In the 
supersymmetry-preserving DR scheme, the Konishi primary is defined as the trace over all 
Ng = 10 — D = 6 + 2e scalars. While the expression (1.4) can easily be modified to sum 
over this regularization-dependent number of scalars, the expression (1.10) is only valid for 
Ng = 6 and hence it is not the highest-weight state of the superconformal multiplet. In 
particular, its anomalous dimension is not given by (1.6), which is usually calculated using 
a descendent of the Konishi operator in the SU(2) or SL(2) sector.*? 

A priori, this discrepancy requires one to abandon the unitarity techniques that employ 
the MN = 4 on-shell superspace. Fortunately, this is not the case. In the following, we will 
show that — at least for the cases at hand — the strictly four-dimensional result can be 
lifted to the D-dimensional result with a simple prescription. 


4.2 From Kz to kK 


Consider a generic multi-loop diagram contributing to the two-point form factor of the 
operator tr(¢;¢7) with outgoing scalar fields ¢x% and ġz. It can only have one of the three 
types of R-charge flow depicted in figure 9 together with the respective tensor structures. 
Only in the case (c) an externally closed scalar index loop exists. The BPS operator 


pı I K I K I K 
p2 J L J L J L 
(a) OrK OL (b) ôrLôJK (c) ôrjôKL 
Figure 9: According to R-charge conservation, only three different contractions of the 
scalar flavors can exist in a generic multi-loop diagram with incoming operator tr(¢r¢,) 


and outgoing scalar fields g and ġz: (a) 67K 671 (blue), (b) drrd7K (green) and (c) 770K L 
(red). 


tr(¢(r¢s)) defined in (1.2) obtains contributions from the cases (a) and (b) but not from 
case (c). The Konishi operator tr(@7¢,) defined in (1.4) obtains contributions from all three 


32The additional 2e scalar components in the D-dimensional continuation of the Konishi operator are an 
example of so-called evanescent operators, which also appear in the context of QCD [78]. See [74] for a 
textbook treatment. 
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cases. The contributions it receives from the cases (a) and (b) are identical to those of the 
BPS operator since the coefficients of the tensor structures in figure 9 do not depend on 
the R-charge. Thus, we can isolate the contribution from case (c) by subtracting the result 
for the BPS operator from the result for the Konishi operator. In D = 4 — 2e dimensions, 
the single externally closed scalar index loop present in the case (c) should not generate 
a factor 6 as found in the FDH scheme but instead the factor Ng = 10 — D = 6 + 2e as 
prescribed in the DR scheme. Hence, in order to obtain the result for K in the DR scheme 
from the one for Kg in the FDH scheme, we simply have to multiply the contributions of 
case (c), i.e. the difference of the Konishi and BPS case, by the ratio 


Nọ 6 + 2e 
=a Te 


(4.1) 


Similar arguments are valid for the three-point form factor of the Konishi operator. In 
our calculation, only its components with either two scalar legs and one gluon leg or two 
fermion legs and one scalar leg appear. While in the former case the previous arguments 
directly apply, in the latter case a slight modification is necessary as shown in the following. 
A generic multi-loop diagram of the latter type, which has incoming operator tr(¢7¢,7) and 
outgoing fields ¢x%, Ya and wp, can only have one of the three possible R-charge flows 
show in figure 10. In addition to the Kronecker 6, the tensor structures contain the flavor 
(or) AB matrices obeying the Clifford algebra. An externally closed scalar index loop exists 


(a) drK (07) AB (b) (or) apo7* (c) dr7(oK) AB 


Figure 10: According to R-charge conservation, only three different contractions of the 
scalar flavors can exist in a generic multi-loop diagram with incoming operator tr(¢;¢,), 
outgoing scalar fields x and outgoing fermion fields %4 and wg: a) 67K (77) AB (blue), b) 
(or)ABOsk (green) and c) 677(0K) AB (red). 


only in the case (c). In analogy to the case of the two-point form factor, we can isolate 
this case by subtracting the result for the BPS operator from the result for the Konishi 
operator. Then, we modify the number of scalars by multiplying this difference by rg. 

The above arguments can be generalized to any number n of points. Moreover, they 
also hold when reexpressing the real scalars ¢; in terms of the complex scalars ¢4pg. This 
allows us to perform the calculations in the FDH scheme, using the superfields (1.7) of 
N = 4 on-shell super space, as done in the previous section. We formulate the modification 
necessary to obtain the correct results of the DR scheme explicitly below. 

In the previous section, we have split the form factor ratios of Ke as 


ae = A nt A , (4.2) 
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where fgpsyn coming from (a) and (b) in figures 9 and 10 is the part identical to the BPS 
form factor and ficgn coming from (c) is the part unique for the operator (1.10). To obtain 
the form factor ratio for K 


fen = ferent Fen (4.3) 
we apply the replacement rule 
om > Téfton =n» (4.4) 


where rg is defined in (4.1). According to our discussion, it should be valid to all loop 
orders.°° 

We have focused on the form factor of the Konishi operator. A similar discussion 
should also be applicable to other operators containing a contraction of flavor or vector 
indices. As in the Konishi case, it is essential to be able to formulate the results in terms of 
two parts, one that contains an externally closed index loop and the other that does not. 
The part without externally closed index loop should be independently computable, such 
as the BPS part in the Konishi form factor. Given such a decomposition, one can then use 
the efficient on-shell techniques, together with a simple modification rule as (4.4). Another 
example of an operator with contracted flavor indices is tr(¢@z7¢r¢K), which has one-loop 
anomalous dimension 8. An example with contracted vector indices is tr(D“¢12D,,¢12), 
which has one-loop anomalous dimension 12. In the latter case, the differences in the 
one-loop two-point form factor between intermediate states in D = 4 and D = 4 — 2e 
dimensions are precisely given by the rational terms in the PV reduction formula (C.5). 
Similarly, i and oo differ by rational terms that are introduced by the replacement 
(4.4); these rational terms arise when multiplying the +-pole from the bubble integral with 
the term in rg that is linear in €. 


4.3 Final Konishi form factors 


Finally, we list the non-vanishing results for the Konishi form factor, which will be used as 

the input in the next section to calculate the cross section. They can be obtained by using 

the form factors computed in section 3 and integral results in appendix B. Note that the 

fk parts have been modified by the rg factor according to the prescription (4.4). The full 

form factors can be obtained as ae = oe + fe A In the following equations, (45) is 
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33 This statement relies on the validity of the DR scheme, which, however, is known to have inconsistencies 
at higher loop orders [79-82]. 
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Two-point two-loop 
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Kslevsa) e(1 — 2€) = + — =q? (s13 + s23)? \—si2 
12 2 € 2 e (4.7) 
a rem aera S) -e(S) I} 
812(813 + 523) (2— = $13 + $93 —$12 —q 
2 2 
z1) eT 812 — $13 (- H ) $12 — $23 ( H ) 
1 a) i a 
I s.6) = e(1 — 2e EE ( ý sig + eo) — $93 ğ 7 S12 + $937 \— 813 


2 
S12 — $13 S12 — $93 MNS 
_3 + se) (4) } — 3r,FB P2,P3,P1,—4) ; 

(cant +513 812 + $23) \—q ia | 


where cr is given in (B.5), q? = s12 + 823 + s31 and we have defined the (rescaled) finite 
part of the one-mass box integral as 


FB(p1, p2, P3, —4) 
-i beet Pe 2 (4.8) 
--8[(4)(-2)+(4-2)-(4)"-24) 


with h(x) = 2Fı (1, -e,1—e,x) — 1. 


5 BPS and Konishi cross sections 


In this section, we compute the cross section discussed in section 2. We first discuss in 
detail the case of the BPS operator (1.2) as a warm-up example. Then, we compute one 
of our main S g oae cross section to two-loop order. We will use the Konishi 
form factors A and o n given in subsection 4.3 that were obtained from the form factors 
of Ke pated in section 3 by applying the prescription of section 4. 


5.1 BPS cross section up to one-loop order 


As a warm-up, we first consider in detail the cross section corresponding to the imaginary 
part of the two-point correlation function of the BPS operator tr(#?,) and its conjugate 
tr($34), (0| tr(?y) (x) tr(d34)(0)|0). Since the operators are protected, the cross section has 
no UV-divergent loop corrections. Moreover, finite corrections do not occur either [83, 84], 
i.e. 

Opps = gha + O (e). (5.1) 


We check this explicitly up to one-loop level. 
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Figure 11: Tree-level squared matrix element. 
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Figure 12: One-loop squared matrix elements. 


Tree level 


Let us start with the tree-level cross section. The squared matrix element, as shown in 
figure 11, is the product of two two-point tree-level BPS form factors, one for tr(¢?) and 
one for its conjugate tr(¢2,). The tree-level non-color-ordered BPS super form factor can 
be obtained from (2.11) and (3.2). It is easy to perform the color factor summation and 
the fermionic integration. This yields the squared matrix element 
2 
Maza D | atm atm ARDANI = =. (5.2) 


a1 ,a2 


The tree-level cross section is given by the integral (E.2) of Mos over the two-particle 


phase space in D = 4 — 2e dimensions. This yields 


2 2 
0) _ 0) (EENS 1 Né—1 
= [APs Moras 7 (7) A(16m)2-*T(Z—e) 2 a 


One loop 


The one-loop cross section is given by the sum of a two-particle and a three-particle channel, 
as shown in figure 12: 


1 
obes = fas: MY, 0 + = J Ps Mes 7 (5.4) 


Two-particle channel 


The squared matrix element of the two-particle channel corresponds to the first two graphs 
of figure 12. As an equation, it reads 


Mis = a > ain dfn Faces Fis 7 Fos P = aM vs R( a , 
aı,a 
™ (5.5) 
where R denotes the real part, and fo. is the ratio between the ¢-loop and tree-level 
n-point form factor of the operator O as defined in (3.1). The tree-level form factor is 
absorbed into Mos For short notation, we denote Fo(1, ..., n) aS Pon- 
There is an important point related to the i0-prescription to be explained here. The 


two-point form factors acquire a factor of (—q? — i0)~¢ for each loop. The function race 


— 35 — 


is the complex conjugate of is and can be obtained from the latter by replacing (—q? — 


i0)~© with (—q? + i0)~€. The sum of both terms amounts to taking the real part of ioe 


Hence, we need the real part of (—q? + i0)~*, which for q? > 0 is given by (see e.g. [85]) 


r +z) z) 


a + 40)? = T(1 + 2x)T(1 — 2a) 


(P). (5.6) 


Using this result to determine the real part of the form factor (4.5) and then inserting 
it into (5.5) together with the tree-level result (5.2) and performing the two-particle phase 
space integral (E.2), we obtain for the first term in (5.4): 


2 
HÎNE 4 Tm 
ee 2 = = f aps, MY) 2 = -= o®, (3) (-3 + T) + O(e) ; (5.7) 


Three-particle channel 


The squared matrix element of the three-particle channel is given by the last graph of figure 
12. The MHV and NMHV tree-level three-point non-color-ordered form factor (2.11) can be 
obtained using (3.2) and (3.7). Performing the color summation and fermionic integration, 
we find the squared matrix element 


0 1 ~MHV,(0) ~*NMHV,(0 ~NMHV,(0) *MHV,(0 
M a= 3] S fan dên dfn Ea fe k + Feps3 ! Ê BPS,3 (a 
a1,a2,43 
2 2\2 
= 2 42,,Ne (IN? — 1) 


= 7g Fr 
aor 812823831 
(5.8) 


Performing the three-particle phase space integral via (E.3), we obtain for the second term 


n (5.4): 
1 BRNE 4 TmT 
aua = = | Ps Ma = of). (=) (3 = T) + O(e). (5.9) 


Summing (5.7) and (5.9) together as prescribed by (5.4), we see that both contributions 
cancel and hence that (5.1) holds at one-loop level. 


5.2 Konishi cross section up to two-loop order 


Next, we compute the Konishi cross section. We start with the discussion of an important 
simplification for the computation, which exploits the fact that a part of the Konishi cross 
section is identical to the BPS cross section (5.1), which is protected. 

At tree level, the squared matrix elements of the Konishi and the BPS cross section 
satisfy the following simple relation*4 


Men= > 2, Piafin = 6D) o Faron Farin =6MBP m+ 5 46) 


colors spins colors spins 
helicities helicities 


34Note that the tree-level Konishi form factors with specified external legs are identical to the correspond- 
ing BPS form factors. 
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where the factor 6 originates from the contribution of all scalar flavor degrees of freedom 
in the two-point function of the Konishi operator (1.4), which does not occur for the BPS 
operator.®° 

Furthermore, the loop correction to the Konishi form factor can be written as linear 
combination of two contributions as defined in (4.3): one that is identical to the BPS 
form factor and the other that is unique for the Konishi operator. We can introduce a 
corresponding squared matrix element that includes a subtraction of the BPS part as 


~ (2 £ e 
Me, = My, E a Ps n ’ (5.11) 
where the factor 6 takes into account that at any loop order £ the contribution Msn 
built from two BPS-type components of the Konishi form factor receives a factor 6 as in 
(5.10). 
Since the BPS cross section (5.1) receives no loop corrections, 


+2 e+1 
oO) =Y ge) f dro = ye f deS AKTS gazi. Gly) 


n=2 n=2 


Note in particular that Me, = 0. Hence, the sum over n can be terminated already at 
+1. As we will see, this simplifies the computation dramatically. 

From (5.10), it immediately follows that the tree-level cross section for the Konishi 
operator also contains an extra factor 6 compared to the one of the BPS-operator, i.e. 


o® = f JPB M, = 6o (5.13) 


At loop-level, it is convenient to factor out o., 


5.2.1 One-loop result 


The bare one-loop Konishi cross section receives contributions from products of tree-level 
and one-loop two-point form factors and of tree-level three-point form factors as shown in 
figure 12. The squared matrix element of the two-particle channel is given by 


1 1 ~(1) (0 ~(0) >(1 0 1 
My’, =J] > ain d'm ( + i 2) = NRU és) , (5.14) 


a1 ,a2 
iati D 0 36 
where we use the abbreviation fk (4 4) = JK oh diay 2ps): 


35The additional prefactor compared to the BPS result comes from the trace over the scalar degrees 
of freedom, which here (in the FDH scheme) is 6. In the DR scheme, one would have to replace 6 by 
No = 6 + 2e. In any case, this factor cancels out when the cross section is divided by the tree-level cross 
section as e.g. in (2.14). 

3The tree-level two-point form factor must contain two external scalar legs to be non-vanishing. There- 
fore, it is not necessary to consider one-loop form factors with other external states such as eae J2- h 
which are moreover zero as shown in section 3. A similar argument applies also to the following two-loop 
computation. 
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Figure 13: One-loop correction from one-loop renormalization constant. 


As discussed above, the result for the three-particle channel cancels with the BPS part 
in the two-particle channel. Therefore, we can subtract the BPS part from (5.14), as in 
(5.11). This yields 


~(1 0 a(1 
Me OM RU is (5.15) 


where sae $6) is given in (4.5). Performing the two-particle phase space integral (E.2), the 


one-loop bare Konishi cross section reads 
m 2\ 6 12 
o% = fas, Me, = of) (=) (-2 — 28) +O(e). (5.16) 


The divergence in (5.16) has to be canceled by the one-loop correction of the Konishi 
operator obtained from the one-loop term Zi) in the operator renormalization constant 
Zg. As shown in figure 13, z% contributes as 


TAES E (5.17) 


Requiring the sum of (5.16) and (5.17) to be finite, we immediately find 


(5.18) 


Comparing this result with the one-loop term of the expansion (2.3) reproduces the known 
one-loop Konishi anomalous dimension yË l= 12, which was first obtained in [45, 46]. 
The renormalized one-loop cross section is hence given by 


2 
En = o% a M = of) (121g 2 = 28) +O(e). (5.19) 


As predicted in (2.15), the coefficient of log also reproduces the correct one-loop anoma- 
lous dimension. 


5.2.2 Two-loop result 


The two-loop cross section is obtained from the contributions to the squared matrix el- 
ements depicted in figure 14. As discussed at the beginning of this subsection, we can 
neglect the contribution that is proportional to the BPS cross section. In particular, it 
is not necessary to consider the contribution from the four-particle channel in figure 14c, 
which involves the complicated four-particle phase space integral. This simplifies the com- 
putation significantly. In the following, we separately compute the contributions from the 
two-particle and three-particle channel. 
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(c) 
Figure 14: The two-loop bare squared matrix element. 


Two-particle channel 


The full contribution of the two-particle channel consists of the three terms 


f smg, + f smg ey + f arsa mg y (5.20) 


where the first term is the bare contribution, and the second and the third term involve 
the one- and two-loop contributions of the renormalization constant Zg, respectively. We 
compute the first two terms. After having considered also the three-particle channel, we 


2) 


will determine zË from the condition that all divergences are canceled. 
The squared matrix element obtained from the bare form factors is shown in figure 


14a. In analogy to (5.15), the first two graphs yield 
~ (2), 0 (2 
M SI R as (5.21) 


where cae is given in (4.6). 

The third graph of figure 14a has no lower-loop counterpart and needs to be discussed 
in detail. It is the product of two one-loop Konishi form factors, and each of them is a 
linear combination of the BPS part ace and the additional i part. After subtracting 
the product of two BPS parts, we obtain 


MET = Me, [r(e ~ TEEN P (5.22) 


where the form factors are given in (4.5).37 

Integrating the sum of the two previous contributions over the two-particle phase space 
(E.2) yields the bare cross section of the two-particle channel with its BPS part subtracted. 
It explicitly reads 


~(2 ~,(2),I ~,(2),II 
5, = f aps. (Ay + AR) 


2\2e[48 184 584-567? 668 
sat(e | he a = INES aA 
K e3 e2 € 3 


(5.23) 


37Recall that for the Konishi form factors we have to specify the external legs while for the universal BPS 


ones this is not necessary. 
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Figure 15: The contribution from the one-loop renormalization constant to the two- 
particle channel. 


Next, we consider the contribution involving the one-loop renormalization constant. It 
is shown in figure 15 and leads to the squared matrix element 


Meg = MED [AR( Aloo) ZO + (ZL) - (5.24) 


Inserting the explicit expressions (4.5) into (5.24) and performing the two-particle phase 
space integration (E.2), we obtain 


2 2 
oe ga = f APS: Mo, 


2\€ 48 144 336 — 287? 
=o|(5) (- 5 - AF -  ~ 70 + gar? + 11265) (5.25) 
36 12 
a H + O(e) 
€ 


Three-particle channel 


There are two contributions to the two-loop cross section in the three-particle channel: 


=| 1 dPS3 MẸ? +f APS Mars (5.26) 


The contribution involving the bare form factor is determined from the diagrams of 
figure 14b. The resulting expression reads 


a ),MHV NMHV , ¢(0).NMHV g«(1—£).MHV 
Mk} = PGD E AS ” K3 Fk 
= (5.27) 


(1) 512 813523 
= 6M Esa) E 7 5 + 2R( A toru,0)) oak 


In the second line, we have not indicated the MHV degree, since the loop correction is the 
ei for the MHV and the NMHV form factor. This allows us to use the abbreviation 


= f9; ),MHV _ p(1),NMHV 


fee = E3 for any fixed three-particle final state. Moreover, we have ab- 


breviated the form factors of the two different final states as fe, re FP (lias Zoss 3g) 

and ‘isha = = f(y pa: 3634) Since the one-loop corrections for different external 

states differ from each other, as given in (4.7), we need to treat the contribution of these 
s23 

two form factors separately. The factors Gh and Tor stem from the squares of the 


(0) 


corresponding tree-level form factors divided by the tree-level matrix element My3- 
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Figure 16: The contribution from the one-loop renormalization constant to the three- 
particle cross section. 


After subtracting the BPS part, we find 


(5 28) 


ge 
(1) (0) ra $12 

Mýs = = 6Mx3 e(z T (2 2)2 + E (p,p, a) (q2)2 
Inserting the explicit results (4.7) and performing the three-particle phase space integration 
(E.3), we find that the contribution to the cross section is given by 

2\2e[ 48 112 224-567? 632 
o a] -5R gy a sasl + 010). 

€ 


(5.29) 
The one-loop renormalization constant (5.18) contributes as shown in figure 16. The 
squared matrix element reads 


MeeVee (5.30) 


Z20K,3 


The corresponding contribution to the cross section can be computed as in (5.9), and it is 


given by 
2 2 
(2) _ 1 (1) _ (0) (LK NE 48 287 
TZ2K,3 = 5 | Ps Maing = OK (=) È — << = 400¢3 + O(e) š (5.31) 


Summing (5.23), (5.25), (5.29) and (5.31), we find 


(2) (2) 
a +o! De + OK 3+ Tz)K3 


_ ,@ | _ 36 of g? 2 
= 0k Ma AT log“ —5 — 384 log ~> + 508 + 727^ + 336¢3| + Ole) . 
e æ H m 
(5.32) 
Since the one-loop UV subdivergences in the bare contributions are canceled by the second 
and fourth terms, all divergences in the above result originate from the two-loop overall 
UV divergence.*® 


Two-loop renormalization constant 


The two-loop overall UV divergence has to be canceled by the third contribution of (5.20), 
which involves the two-loop renormalization constant. In analogy to (5.17), this contribu- 
tion reads 


= = 20.0) 2) (5.33) 


38 All infrared divergences should be already canceled between the different channels. 
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Requiring the sum of (5.32) and (5.33) to be finite, we immediately find 


@ 18 12 
Ln == 


= (5.34) 


€ 


Comparing this with the expansion of (2.3) in terms of the anomalous dimension to two- 
loop order yields the known one- and two-loop Konishi anomalous dimension first obtained 
in [45, 46]: 

P =12, 7 = 48. (5.35) 


Final result 


Summing (5.33) and (5.32) yields the renormalized two-loop cross section 
D _ (0) q q? 
Trh ET 72 log? =, — 384 log = + 508 + 727° + 336¢3| + O() . (5.36) 
' H H 


Finally, we compute the second order term in the expansion of (2.15) and obtain 


TK,R Á or 1 EER 2 q? 2 
log ( R = AR __( G8) = —48log = +116 + 720? + 336¢3+O(€) . (5.37) 
KO 50) 2K 40) u? 
K K K 


We find that the coefficient of log 5 gives the correct two-loop anomalous dimension, as 
expected from (2.15). 

Including also the one-loop result (5.19), the logarithm of the normalized Konishi cross 
section is given by 


2 2 
log (=) =g? (2 log 2 — 28) +g! (-48 log Z +116 +727? + 336s) +0(gf, €). 
Oo 


(5.38) 
The finite terms that are independent of log s yield the constant C, and by a comparison 
with (2.16) they determine the one- and two-loop terms of the constant M in (2.4). In 
particular, this yields the full two-point function (2.4) up to two-loop order. 


Some discussion 


The renormalized cross section can be computed in different ways. In the above presen- 
tation, we have treated the bare contribution and the terms involving the renormalization 
constant separately at the cross section level. One can also first perform the renormaliza- 
tion of the operators via form factors, as described in appendix F, and then compute the 
renormalized cross section directly from them. Furthermore, the terms involving the renor- 
malization constant can be obtained directly by expanding relation (G.1), which gives the 
renormalized cross section in terms of the bare one and the renormalization constant. For 
example, the sum of (5.25) and (5.31) that involve the one-loop renormalization constant 
can be obtained as?’ 


2 2 1 1 1 0 
Giga besiege oe FEL) ce (5.39) 


3°Tn this case, one needs the result of the one-loop cross section up to O(e) order. The result (5.16) based 
on (5.12) is not enough, since opps is not zero at O(e) order. 
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We have checked that these different ways give the same result. 

As discussed in appendix G, the above result depends on the renormalization scheme. 
One can define the new coupling at which the subtraction is performed as gg = ge% 
and then expand the expressions in terms of the original coupling g. This scheme change 


e?el: in all the above equations. With such a 


can be implemented by replacing zit = zę ) 
modification in the above computation, one finds that the renormalized cross section (5.38) 
acquires a finite additive contribution 2yo, demonstrating that M in (2.4) is scheme- 
dependent. This agrees with the expectation from (G.7), since the scheme change can be 
understood as a change of the ’t Hooft mass y —> we. 

Finally, let us briefly comment on the FDH scheme we have chosen in the computation. 
In the FDH scheme, we set the number of external scalars to 6 and use polarization vectors 
in D = 4 dimensions for the form factors. As discussed in section 2, this corresponds to 
the prescription given in (2.19), where the sum of the degrees of freedom for the external 
legs is performed by the ņ-integration based on the SU (4) representation. One can also 
perform a detailed analysis at the diagrammatic level, as for the form factors in section 
4, which leads to an alternative prescription for obtaining the cross section in D = 4 — 2e 
dimensions. We will not present the details in the paper, but we have checked that both 
prescriptions give identical results at least up to the two-loop order. 


6 Conclusion and outlook 


In this paper, we have studied form factors of non-protected operators in MN = 4 SYM the- 
ory, specifically of the Konishi operator, using on-shell unitarity techniques. Importantly, 
we have found that this requires an extension of these techniques. We have obtained ex- 
plicit new results of the three-point form factor at one-loop and two-point form factors up 
to two-loop order, given in (4.5)-(4.7). The application of on-shell methods to determine 
such form factors, which are partial off-shell quantities involving both, composite operators 
and on-shell states, provides a step to deepen our understanding of the connection between 
modern on-shell techniques and the off-shell world of correlation functions. 

Another important aspect of this paper is to provide a physical observable within 
N = 4 SYM theory, given by a cross-section-type quantity: the inclusive decay rate of 
a state, described by a composite operator carrying timelike off-shell momentum q, into 
any final on-shell multi-particle state. We gave a formulation of how to compute this 
observable. Using the Konishi form factor results mentioned above, we performed an 
explicit computation of the total cross section up to two-loop order, given in (5.38). Via 
the optical theorem, this also yields the two-point function up to this order. 

The UV divergences appearing in the Konishi form factors together with the IR diver- 
gences require the renormalization of the operator. This is carried out explicitly in the com- 
putation of the total cross section in which the IR divergences cancel. We have reproduced 
the known Konishi anomalous dimension up to two-loop order from the renormalization 
constant and also identified it as the coefficient of the log 5 term in the renormalized cross 
section (5.38). 
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Since the Konishi operator is not protected by supersymmetry, interesting subtleties 
and new features appeared, which we now summarize. 

First, an important subtlety occurs in the unitarity-based computation of the Konishi 
form factors. The Konishi primary is a trace of all scalars. In order to preserve supersym- 
metry, the Konishi operator has to be continued to D = 4— 2e dimensions, i.e. the number 
of scalar field flavors that is summed over has to be Ng = 10 — D. However, working 
with four-dimensional unitarity based on Nair’s on-shell superspace, one can only directly 
compute the form factor for the different operator Kg, which in D = 4 — 2e dimensions 
has Ny = 6 scalars rather than Ng = 6 + 2e. In order to find the Konishi form factor, 
the results based on four-dimensional unitarity have to be modified when they are lifted 
to D = 4 — 2e dimensions where the occurring divergences are regularized. We provide a 
rigorous prescription (4.4) for this modification that yields the form factors of the Konishi 
operator. 

Second, the Konishi form factors contain some interesting features that are not present 
in other on-shell quantities studied in M = 4 SYM theory studied so far, such as scatter- 
ing amplitudes or the BPS form factors as partial off-shell quantities. The bare Konishi 
form factor contains bubble integrals and bubble subintegrals, which are UV divergent. 
Moreover, the one-loop three-point result contains a rational term. In addition, the co- 
efficients of the individual integrals occurring in the form factor results involve spurious 
poles, which only disappear after multiplication with the integrals and summation over all 
contributions. Last but not least, the loop corrections of the Konishi form factors with 
different external states turn out to have quite different structures, even if they are in the 
same MHV sector. The emergence of these features that are familiar from QCD can be 
traced back to the fact that a non-protected operator has been inserted into the action, 
formally breaking its supersymmetry. 

Finally, let us briefly mention some further directions one can pursue following this 
work. 

First, it should be straightforward to generalize the computation of the one-loop Kon- 
ishi form factors to the higher-point cases. It is also interesting to proceed to higher loop 
orders. In particular, using the known IR exponentiation property of the Sudakov form 
factor, the knowledge of the two-point Konishi form factor alone is sufficient to extract the 
Konishi anomalous dimension. We explain this in appendix F. Turning the logic around, 
we also give a prediction for the three-loop two-point Konishi form factor apart from finite 
terms there, only using the known three-loop anomalous dimension and the IR exponenti- 
ation in addition to our form factor results. 

Second, our detailed example of how to apply four-dimensional unitarity to compute 
the Konishi form factor by understanding an encountered subtlety and providing a solution 
is a solid stepping stone for further studies of other non-protected operators, based on 
generalizing the prescription we give in section 4. Combining our insights with those from 
the recent one-loop calculation in [33], it should be possible to compute the minimal form 
factors for general operators at two-loop order via on-shell methods. This would allow us 


404 similar rational term also occurs for the minimal form factor of operators in the SL(2) subsector [33]. 
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to determine the complete two-loop dilatation operator of N = 4 SYM theory which yields 
all two-loop anomalous dimensions as eigenvalues. Besides the anomalous dimensions, the 
other important CFT data is given by the structure constants, which can be computed 
from three-point functions. It would be very interesting to use similar unitarity-based 
techniques to compute them. 


Furthermore, as given in [25, 38-40], the so-called energy energy correlation function 
can be computed as a weighted cross section which is very similar to the total cross section 
studied in this paper. In [39, 40], different techniques not relying on cross sections have 
been used to evaluate them up to two-loop order. It would be interesting to obtain them 
also from direct cross section computations. The interpretation of the cross-section-type 
quantities at strong coupling via the AdS/CFT correspondence is also an open problem. 
In particular, it would be interesting to consider the phase-space integration with strong 
coupling form factors in the framework of string theory. 


Finally, as a cousin of N = 4 SYM, the so-called ABJM theory [86] has been intensively 
studied in recent years. In particular, the form factors for half-BPS operators have been 
determined in this theory as well [87-90]. It would be interesting to pursue a similar study 
as in this paper for the ABJM theory, especially for the form factors of non-protected 
operators. 
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A Fourier transformation of the two-point function 


In this appendix, we give our conventions for the Fourier transformation and for the tran- 
sition from Euclidean to Minkowski signature. 

In Euclidean signature, the coordinate dependence of the two-point function has the 
following representation in terms of the momentum space integral: 


1 = D248 (4 = | dP qe eE TE 
(z3)^ r(A) (27)? (q2)7-4 


(A.1) 


D-1 
where qp : £p = qe, 0%, 0 + Da qiti. 
In Minkowski signature with a mostly-minus metric, the exponent in the momentum 


space integral is given by —iq - x, and the integrand has poles at gg = +\g¢|. We want 
positive energies gg > 0 to propagate into the future xp > 0. Hence, the pole at q = |q| 
has to be picked when for xp > 0 the integral over qo is closed in the negative imaginary 
half plane such that the exponential factor vanishes for qq + —ioo. This is guaranteed if 
we replace q2 — —q? — i0 in the denominator of the above expression. The position of the 
poles fixes the Wick-rotation to be counterclockwise in momentum space, i.e. go = iqr,0, 
and for x9 = —i%p,9 to be clockwise in configuration space. This leaves the exponential 
invariant, and it can be transformed to Minkowski signature by flipping the sign of the 
spatial momenta qi. We hence find 


1 B Dr(2 _ A) dPq ena 
arrays OP ME faa gigs AF 


D- 
where q: £ = qozo — De qiti. 
B Feynman integrals 


In this appendix, we present all integrals that enter the form factor results in section 3, as 
well as our conventions. Moreover, we show how the cut integrals are lifted to full integrals. 

As a regularization procedure, the four-dimensional M = 4 SYM theory can be con- 
tinued to D = 4 — 2e dimensions. Both IR and UV divergences are then captured as poles 
in e. Moreover, the Yang-Mills coupling constant gym has to be replaced by gymp, where 
u is the ’t Hooft mass which is introduced in order to keep gym dimensionless [91]. Hence, 
the combination gu‘ with g given in (1.5) occurs as the effective loop expansion parameter 
of the perturbation series in the large-N. expansion. 

From Feynman diagrams, the following combination of the coupling constant, the ’t 
Hooft mass u and the loop integral occurs at -loop order 


ont ay f Aah dlp f(li,..., le) 
(9ymbt ) Ne (—i) f (27)? aa G La- 


where the integral J is of the following form: 


= g1O[f(ly,...,le)] , (B.1) 


dP? ly dP ly fl, as , Ip) 


zai: rs 27 pede 
a7 2 am 2 I]; D; 


1OLf(I1,..-5le)] = (e™ p)“ f 
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In these formulae, the D;’s are the propagators, i.e. Dj = k? +120 for k; being the combi- 
nation of external momenta and loop momenta that flows through the propagator. 


Lifting the cut integral 


Let us explain our procedure and conventions for lifting the cut integrals to the full integrals. 
Consider the triangle term in (3.15) as an explicit example. We have 


l pı pı 
2 $12 2 ; -2 
Ne | dPS —_—— = Nes — —ig’s ; 
Ivm a 20} (h T pı)? Same voaa s p lifting g S =< b 


l2 
(B.3) 
where the phase-space integration measure dPS» 17) is defined according to (2.20) with 
dPl š š x 
measure factor TaT: To lift the cut integral to the full integral, two steps have to be 
performed. 


One is to replace the cut propagator as 


i 
2784 ($) — z (B.4) 
j 
The other is to change the measure factor and coupling constant as in (B.1), such that 
the integrals of uncut graphs are defined in terms of (B.2). This prescription is employed 
throughout section 3. 


List of integrals 


In order to present the expressions for the required integrals in a compact form, we define 
q = J; pi and introduce the factor 


r-T +e) 


= e%E€ 
° Tl — 2e) 


cr (B.5) 


Moreover, in the results for the integrals, all (—q?)“ should be understood as (—q? — i0)” 
and similarly for (—s;;)“. 


In the conventions introduced in (B.2), the one-loop integrals that are required to 
compute the one-loop form factors read 


dl 1 ao jay 
r pS a E B. 
X sey o tae GF) a 
pı 
dP] 1 a 1 -fj 
= YE ,,2\€ e E S B.7 
[_ (e u°) |= (1 + p1)? (l — po)? e2 ay ( a n ABO 
Pı D 
<T = (e™ p’): a. 
a in? (l+ pi t+ p2)? (l — ps) 


ar A ~ (FF) 
= —- > — | — B.8 
E€ 813 + S23 pe p? ee) 


—AT— 


. q?1 1 
ps = ee 0 f Spa or B.9 
<S = (ey | PPE +p tn +m @ vr 


02N E 
= cr 2 [ = (=) oF (1, zei E= qĉs13 ) 


€ 819893 u? 812823 
— 812 —eE 
+ ( ye ) 2Fı (1, —€, 1 — €, —) 
— S23 —eE 
+ ( ie ) 2Fi(1,—e,1 — e, —22)| , 


where 2F} denote the Gaussian hypergeometric function. The results (B.6)-(B.9) can be 
found for example in [92-94]. 
Furthermore, we need the following integral, which evaluates to a rational term [92]: 


peo em ity fEl n 
= (Q —_ E 
ODE ETETE ih 


oe ae 
(1 — 2e)(2 — 2€) 513 + 593 u? u? ? 


where the momentum le in the 2e dimensions is introduced in (C.1). 
To compute the two-loop two-point form factor, we need the following two-loop inte- 
grals. Using IBP identities as e.g. implemented in LiteRed [95], they can be reduced to 


master integrals ast! 


3(1 — 2e)(1 — 3e PY (1 — 26)? ae 
ote Pac +! H(X S). 
l Pı Pı 
© (2 — 3e)(2 — 9e + 10e? — 4e?) 
$1152] -<C s = 0-A- AeA <x r (B.13) 
1— 4e + 2e o 2— 364+ 2 a 
“ee “a 2(1 —e)e (<< ) , 
l pı pı 
_ (1 — 2e)(2 — 36)(3 — 5e) 
$1152] << : = E —4e(—-@) LX 7 (B.14) 
= (1+.e)(1 — 26) „N 7 e(—q?)? << 
e(1 — 4e) p (1-46) p2 í 


“Recall that the loop-momentum-dependent prefactors are understood to appear in the numerators of 


the depicted integrals. 
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where the master integrals are [96] 


ed T(1 — e)?T(1 + 2e) q?\~** 
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3F2 (1,1,1 +26; 2 + €,2 + 26; 1) 


3P (1, 4e, —2e; 1 — 3e, 1 2e; 1) 


-H 


4F3(1,1 — €, —4e, —2€; 1 — 3€, 1 — 2e, 1 21)| 


(B.15) 
and the one-loop bubble integral is given in (B.6). 


C Passarino-Veltman reductions 


In this appendix, we summarize some results on Passarino-Veltman (PV) reduction [97], 
which we need in section 3. 

We use the four-dimensional-helicity (FDH) scheme of [76, 77], and decompose the 
D-dimensional loop momentum / into a four-dimensional part la) and a (D — 4) = —2e 
dimensional part le, where we assume that e < 0. This yields the following decomposition 
of the scalar product: 

qal = Vy =P +E, (C.1) 


where Nuy is the four-dimensional metric.*? Arbitrary four-dimensional external reference 
momenta are denoted as ki. 

Bubble. The D-dimensional bubble integral with external momentum q defined in (B.6) 
may include a non-trivial polynomial f(l) of the loop momentum / and the reference mo- 
menta k; in its numerator. Denoting this integral as IP [f (D](q°), we find the following 
relations for the reduction of tensor integrals to scalar integrals: 


PEND = -EP Pe), (C2) 
IP kE kay?) = (LANE) _ EO) ppg 
Z (£ ` — ` k2) -_ a) ie lia") ; (C.3) 


“Note that, although D = 4 — 2e, there is a plus sign in front of l? since the metric has mostly-minus 
signature. 
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Triangle. Next, we consider the D-dimensional triangle integral with numerator f(l), 
which depends on two arbitrary momenta qı and q2. It is defined as 


d?I fÙ 


PFO (a, ge [5 ei C4 
x [F], q2) = (e® u^): ime PUF gl q)? (C.4) 
We find 
2 D, 2 aig 
IPK kDa) = > S 5 2 Ra , 
i=0 i=1 (C.5) 
IP [( - k1)(l - k2)\(q1, 42) a © + C3o1P + C3 dp [i |(q1,42) , 
where 12°) = IP (u ~ 2)"), I = 1P (8), Ip = 1P), 
TE i 2 
C= -7(So ai @) azı — 7011 (42: k2) + (k k2) — 0] 1 2 (A Jai, 
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Of =—2( aug?) a22 — taalan ke) + Elk ka) — 8)42 (AY 
: AN k a ° ! i=l k 
(0 a, 1 (C.6) 
Cy’ =-) 0; + z(k1: ko) , 
2 2 i 
2 2 2 21 
C30 = (Zud) (Zeug) + li (ki- ke) ap aA, 
C3 e = Mieli 
with 3 
m=1 
2 
ag = — Xai ; Qi = ai; (C.7) 
i=l 
2 a: 
b= Do (aA ulg k), A= > Ay 
j=l j=l 


The integrals involving l? give rational terms [92]. The rational term for the rank-two 
tensor triangle integral is given in (B.10). 


D Checks of the three-point one-loop Konishi form factor 


Rational term in FO (14,95 2434) 3g+) 


An interesting feature of the Konishi form factor is the occurrence of rational terms at one 
loop. 
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Figure 17: Feynman diagrams of the one-loop three-point Konishi form factor that con- 
tribute to the rational term. 


For the form factor FO (1419s 2dsar 8 34+), denoted as Fy: ( this corresponds to the 


) 
K,($,6.9)’ 
triangle integral containing the /?-term; see (3.45). Using (B.10 


0), the rational term, denoted 
by R{e], can be computed as 


j= Fo) No 813823 


(1) 
F — mM. 
it K,($:99) s19(813 + 523) 


K,(:6.9) (D.1) 


Since the computation in section 3 is based on the four-dimensional unitarity method, one 
might be concerned whether additional rational terms are missed in this approach. In the 
following, we show that the above result is actually complete by comparing with a Feynman 
diagram computation following the strategy of [98]. 

First, from the power counting criterion given in [4], a one-loop integral can generate 
rational-term contributions only if it has a high enough power of the loop momentum / in 
the numerator of the loop integrand, which is given by 


m>n-2 for IP{(i)™| withn>2, 


m>1 for IPD"). (D-2) 


Second, we can safely neglect Feynman diagrams that appear in the computation of the 
BPS form factor, since the sum of them is known to be free of rational terms [20]. From 
these findings, it turns out that only two diagrams need to be considered, which are shown 
in figure 17. 

Using standard color-ordered Feynman rules (see e.g. [67]), these two graphs give 


(0) ANg \ (23)(31) ye aye f AP (L- IL-1 — p2) 
FK E “aay O p’) [2 A ns 2(1 + p12)? , aa 
V26X3 


where the polarization vector is given by e = = and € is an arbitrary reference spinor. 


Then, applying the identity” 


PU (1 ki) ko) O kiska (ki - ka) (ko a1) — 
R| f Santer = a nn ie =v; (D.4) 


2 2q1 ` q2 
we can extract the rational term of (D.3) immediately, which, after some simple spinor 


algebra, turns out to be identical to that given in (D.1). Thus, we have shown that the 
unitarity method gives the complete rational terms. 


43This identity can be obtained using PV reduction and (B.10). Formulae for more general cases can be 
found in [98]. 
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Spurious poles 


The coments of the integrals in (3.45) and (3.46) contain unphysical poles, such as the 
pole EEA = ee Such poles cannot originate from propagators in the underlying 
Feynman diagrams. Physical consistency requires that they must cancel when multiplying 
the coefficients with the respective integrals and summing ai contributions.*+ Here, we 


check that this is indeed the case. We focus on the pole the other poles can be 


ET 
treated in a similar way. 

Let us first consider the case of Fiy; 2y,3¢). Only the coefficients of the bubble 
integrals contain spurious poles. Summing over all bubble integrals, the + term is free of 
the pole, and at finite order we find 


- 32i log (22), (D.5) 
S12 + $13 S23 


This is indeed finite for s12 + s13 —> 0, as can be seen from the expansion 


log(1 + x) zc x 3 
al ee ee ee D.6 
- 4E toa’) (D6) 
with z = 124513 


The FO (14,24, 39) case is a little more complicated. In this case, both the bubble 


and triangle integral contain the pole en in their coefficients. Expanding to finite order 
and extracting the terms that contain this pole, we find 
$138 s $138 
_ 13823 ; ( =) 13823 (D.7) 
(s13 + s23) S12 812(813 + 523) 


where the first term stems from the sum of bubble integrals and the second term is the 
rational term. Each term itself is divergent when taking the limit s13 + s23 — 0; however, 
the sum of the two terms is finite in this limit. 


E Phase-space parametrization 


In this appendix, we provide formulae for the parametrization of the phase-space integrals. 
Furthermore, we give details on the non-trivial three-particle phase-space integration en- 
countered in section 5. 

The n-particle phase-space integral is defined as 


I On 54 (o}) (en) PP (Xeo (E.1) 


n 


faso- (T 


where (e) denotes the integrand, i.e. the squared matrix element. 
When n = 2, the squared matrix element depends only on q?, and we can evaluate the 

two-particle phase-space integral independently: 

(°) 


fas; (e) = fps2 (0) , fes2 = 16r)? TG- 6) i 


(E.2) 


“4This is a common feature for one-loop QCD amplitudes, see e.g. [5]. 
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The three-particle phase space can be parametrized as 


1 1 
fes (e) = fes f daa "(1 ae dy [y(1 — y) (e), (E.3) 


with 


(qq 


fess = 2(4r)3 -T (2 — 26) ` (E) 


The ratios of Mandelstam variables occurring in the squared matrix element are parame- 
trized as 
ži Sjk 1 


in which (i, j, k) can be any permutation of (1, 2,3), since the phase-space measure is totally 
symmetric in p1, p2, p3- 


Some details about the three-particle phase-space integral 


The phase-space integration becomes non-trivial for the squared matrix element involving 
the three-point one-loop form factor. It contains the finite part of the box integral (4.8), 
which involves the hypergeometric functions 2f}. 

The corresponding phase-space integrals, which are necessary to evaluate (5.26), are 


1 2 el 
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F Anomalous dimensions via two-point form factors 


In the main part of this paper, we have determined the anomalous dimension of the Konishi 


operator from its cross section, i.e. from the imaginary part of its two-point function. It 
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is also possible to determine the anomalous dimension of the Konishi operator from its 
two-point form factor alone. As seen throughout this paper, form factors of non-protected 
operators contain both UV and IR divergences. To extract the UV divergences, one needs 
to subtract the IR divergences. The computation of the cross section, as done in section 5, 
is one of the safest ways to do so. On the other hand, the IR divergences, in particular for 
Sudakov form factors, have an universal structure [16-18]. This allows us to subtract the 
IR divergences directly from the form factors. The remaining divergences are purely UV 
divergences, from which we can read off the anomalous dimension of the operator.*° 

In terms of the effective planar coupling constant (1.5), the logarithm of any (renor- 
malized and diagonally renormalizing) Sudakov form factor in MN = 4 SYM theory has the 
following structure; see e.g. [19]:*° 


OO OO 2 le 
l = 22 o) G 24 (£) — Ias Og Pe) VO 
og fo.R 2 g (log for) 2 g E le)? Dee + Fin’ ) +0(6), 
(F.1) 
where the pole terms originate from the IR divergences and are determined by the universal 


cusp and collinear anomalous dimensions 


Yeusp(g = Yost = = 89° = 16629" F 176¢49° F O(g) , 
(F.2) 
5 
= SD gPa = —ACag! + 16(265 + = Cal) 9° + O(g”) . 
£=1 


The finite terms of the logarithm of the form factor depend on the specific properties of 
the form factor such as the choice of the operator. In particular, they contain a remainder 
function, which was studied in [27, 29]. 

For non-protected operators, renormalization is required. The renormalized form factor 


is given by 
L 
L-8) 
a n ice (F.3) 
where the renormalization constant is related to the anomalous dimension as shown in 
(2.3): 


The universal structure of IR, divergences, together with the bare Konishi form factor, 
allow us to determine the renormalization constant and therefore the anomalous dimension. 
In the following, we employ the two-loop Konishi form factor to reproduce the Konishi 
anomalous dimension (1.6) up to two-loop order. Reversing the logic, we then give a 
prediction for the bare three-loop two-point Konishi form factor up to and including O(e7!) 
order by using the known three-loop anomalous dimension. 


45-This route was also taken in the unpublished notes of Boucher-Veronneau, Dixon, and Pennington [71]. 

46This form was checked for the minimal form factor of the BPS operator tr(¢i2) up to the third loop 
order [14], for the minimal form factor of the BPS operator tr(¢f2) up to the second loop order [29], for the 
n-point MHV form factor of the BPS operator tr(¢72) up to the first loop order [20] and for the 3-point 
MHV form factor of the BPS operator tr(¢?2) up to the second loop order [27]. 
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One-loop form factor 


The one-loop bare form factor is given in (4.5). From the universal IR structure, we know 
that 


2 (1) (1) 
H E Ycus G 
(log fcr) = Sek = feat Ze = E ~ Te 7% ) + OC), ae 


where the one-loop cusp and collinear anomalous dimensions are given in (F.2). The simple 
pole in i, has to be canceled by the one-loop term in the operator renormalization 
constant, which yields 
1 
-T (F.5) 


in agreement with (5.18) and the known one-loop anomalous dimension. Thus, the one-loop 
renormalized form factor is 
(1) pw? \e 21 +e+e)eIT(—-c)*T(1+e.) 6 
Ikk- ( ) —— r a aA (F.6) 
i (—1 + 26) (1 — 2e) € 

Two-loop form factor 
The two-loop bare Konishi form factor is given in (4.6). From the universal IR structure, 
we know that 


pera = e484) = (1822 +22) -3 (0) 
mp al) (F.7) 
= (S)"(-TE- B) +00. 


where the two-loop cusp and collinear anomalous dimensions are given in (F.2). This 
determines the two-loop term of the renormalization constant as 
Zz -2-2 F.8 
e (F.8) 
which perfectly agrees with (5.34) and yields the known two-loop anomalous dimension. 
Hence, the two-loop renormalized form factor is 


2x2 I28-= 56-r — 2s 267? Trt 
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(F.9) 
Prediction for the three-loop bare Konishi form factor 
Now, we reverse the logic. From the universal IR structure, we know that 
(1) \8 
(log fer)® = fE} - Heath +5 (AEk) 
3) 2) (1 3 2) (1 1/0) \3 
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where the three-loop cusp and collinear anomalous dimensions are given in (F.2). Using 
the known one- and two-loop form factors, and together with the renormalization constant 
up to three loops obtained from (2.3), (1.6) and given by 


22) _ ao ae) le Ik) _ 36 72, 56 (F.11) 
K 4863 8e2 6e e e E” i 


we can predict the three-loop bare Konishi form factor as: 


T? z4 


—q? 36 A e3 e2 
4744 — 141r? — 554¢3 — Sint 4 BTC | 87865 
= G3 40 54 B |4 O( e?) , 
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This should be compared with a direct computation. 


G Renormalization-scheme transformations 


In this appendix, we review transformations between different mass-independent renormal- 
ization schemes and derive the behavior of the cross section (2.15) under such transforma- 
tions. 

A renormalization scheme specifies a regularization procedure for the UV divergences 
encountered in perturbation theory beyond tree-level and a prescription for the subtraction 
of these divergences into renormalized fields, coupling constants and composite operators. 
The subtraction prescription specifies how the UV divergences are removed from the pertur- 
bation series. In particular, it has to be indicated which finite parts are absorbed together 
with the UV divergences into the counter terms or — equivalently — the renormalization 
constants determining the relations between the bare and renormalized quantities. 

A modified renormalization scheme, which contains a different prescription for sub- 
tracting the UV divergences from the perturbation series in g, can be described by applying 
the subtraction of the original scheme but to the perturbation series in a modified coupling 
constant gp = ge®®. Thereby, the parameter o specifies the finite terms that are subtracted 
together with the UV divergences. Since the combination gu‘ of the coupling constant 
g and ’t Hooft mass u is the expansion parameter of the perturbation series, the change 
between schemes, i.e. between g and gg, can easily be implemented by changing u. If we 
demand gou% = gu‘, the transformation of the perturbation series to the scheme ọ, but 
written in terms of the original coupling constant g, is given by replacing u > Ho = pe 8. 

A widely used renormalization scheme is the dimensional reduction (DR) scheme, which 
combines regularization by dimensional reduction with a minimal subtraction of the diver- 
gences into counter terms or — equivalently — renormalization constants. Minimal sub- 
traction means that no finite terms are subtracted. In the DR scheme, minimal subtraction 
is applied to the perturbation series in the coupling constant YA, A= gu Ne- 

In this paper, we work in the modified dimensional reduction (DR) scheme. It employs 
dimensional reduction as a regularization procedure, but the subtraction is non-minimal 
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coupling constant g defined in (1.5). Hence, the subtraction procedures of the DR and DR 


in terms of the coupling constant = g2uNe. It is, however, minimal in terms of the 
scheme are related in the same way as those of the famous MS and MS schemes defined 
in [91] and [44], respectively, that employ dimensional regularization as a regularization 
procedure. The expressions in the former schemes are obtained from the ones in the latter 
schemes by replacing u — pe 2, where o = 3 (log AT — W). 

Consider the renormalized cross section og. Inserting (2.1) into (2.8), it is given as 
the product of the squared operator renormalization constant Zo introduced in (2.1) and 
the bare cross section op: 


oR = Zo(g,€)7op . (G.1) 


The logarithm of the ratio of the bare and the tree-level cross section ¢ then has the 
following expansion up to two-loop order: 


( 


1) 2.9 )y2 
OB _ af Yo. , )\, aN O Yo. , .2)_ (89) 6 


The finite terms 5) become the coefficients of the perturbative expansion of the ratio of 
the renormalized and the tree-level cross section 


2 
aol = aa [1+ gs) + g 8) SOG", €)] (G.3) 
in the DR scheme where the coupling constant is (1.5). 

The condition gou% = gu“ implies that the expression (G.2) is the same in all schemes. 
However, only the subtraction prescription of the DR scheme leads to the expression (G.3) 
for the renormalized cross section. 

The renormalization constant Zo(g,¢) of the DR scheme obtained by performing mini- 
mal subtraction at the coupling constant g can be expressed as the renormalization constant 
Zoo = Zo(go,€) in the scheme o obtained by performing minimal subtraction at the cou- 
pling constant gy times a factor without poles in e. Hence, the difference of the logarithms 
of these constants is finite and given by 


log Zo(g, €) = log Zo,9 + AZ + g6AZG? + O(g) , (G.4) 
where 
g 2, 
AZg® = (() -— 1) Z% = -P o(1 - ge) + O(€) , 
(2) g\4 a (AZ — ay, 2a) (2) a 
AZ% = (5) — 1) 25 — — Zo AZo = -yö 0 +0(e). 
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We have used the expansion given in (2.3). In the expression Az we have kept the term 


) 


linear in €, since it leads to a finite term in e€ in the expression Az’ , when it is multiplied 


by zy). 
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Adding (G.2) and (G.5) leads to the following relation of the renormalized cross sections 
in both schemes 


OR 
log —~ = 


o = lg Ty AP AZy) + G'AZG?) + O(4°,€) = log T+ 2700+ O(9°, e) , (G6) 


gl) 
where we have inserted gg = ge® and neglected terms that vanish when e —> 0. This 
relation can be interpreted in two ways, as follows. 

First, one can insert the expansion (G.3) for op and the same expression for CR,oọ but 
with u replaced by fy. Then, one obtains the relation jug = we ° mentioned already at 
the beginning of this appendix. This shows that a scheme change can be performed by 
changing ju. 

Second, one can insert the expansion (G.3) for og and a similar expression for oR,, but 
with the finite expansion coefficients 5!) replaced by a. Then, one obtains the behavior 
of the finite terms under a scheme change, given by the relations 


a. = 3) + 29 0 : (G.7) 


H Feynman diagrams 


In this appendix, we compute the unrenormalized form factors of section 3 to two-loop 
order via Feynman diagrams. See e.g. [99] for the Feynman rules of the MN = 4 SYM theory 
in our conventions. In particular, we demonstrate how the analysis of section 4 works for 
the concrete diagrams and that we did not miss any rational terms in section 3. 


One-loop self energies 


For the calculation of the unrenormalized two-loop form factors, we need the one-loop 
self-energies of the gauge and scalar fields. They occur as subdiagrams in certain two-loop 
diagrams. 

The one-loop self-energy of the gauge field is determined from diagrams in which the 
scalar fields, the fermion fields, the gauge field itself or the ghost field propagates in the 
loop. They evaluate to 


2 os 
g , a 
COS = FNO” Is wy ’ ate f= g Ny" Lewy y 


(H.1) 
g cos i _ 
A = z’ (DI; w + Iph uv + DT i uv) ’ ey a. =g ô Ion pw , 


where g is the coupling in the DR scheme defined in (1.5), and besides the number of 
scalar flavors Ny = 6 + 2e we have also introduced the number of fermion flavors Ny = 4 of 
N = 4 SYM theory. Moreover, we have split the contribution from the gauge loop into the 
tensor integrals Is, Ign pv occurring in case of the scalar- and ghost-loop contribution, 
respectively, and into Jpn, , which is associated with the remaining physical degrees of 
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freedom of the gauge-field polarizations. The occurring integrals are expressed in terms of 
the simple bubble integral in (B.6) as 


e f Pl (q— 2) u(g— 2), 1 
La = (0 PY f= pap tate — Pe) Oe 


m2 


dl tré,(q—Do,l D-2 
fv = (e Cae bee: Pq? Dy ee T Iw) , 


< faPl 1 
Tyne = EEP) | pag rp tna + taula- Dv + Md = Duar = 8a- Dla = De) 


inZ 
+ (50° — 2q - (q — 1) + 2(q-1)”) gv) 
4D -2 
E Os, 


D — 
Ign wv = (ce wy f : agit _ woe = 2) Guy = guv) =< = . 


m2 


(H.2) 
Inserting the results for the tensor integrals into (H.1) and summing all contributions, 
we obtain 


2 
e- = Fo" (No + D)Is w + 2Nyle yy + Ih wv) 


2 
9? a No + D+ 2(D —2)Ny — (4D — 2) 


= 2976 (ququ — ggv) =~ = 


The first line shows that our decomposition of the gauge-loop contribution in (H.1) is 


advantageous: D and Ng only appear in the combination D + Ng which is insensitive 

to the simultaneous continuation of D and Nọ as prescribed by the DR scheme, cf. the 

discussion in section 4. We note that when inserting the appropriate numbers flavors in 

the second line, the dependence on D originating from the tensor integrals is also canceled. 
The remaining one-loop self energies for the scalar and fermion fields read 


(H.4) 


One-loop form factors 


In the Feynman-diagram approach, the one-loop form factors for the BPS operator (1.2) and 
the Konishi operator (1.4) are obtained from the two diagrams given in table 1. Completing 
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Table 1: Diagrams for the unrenormalized one-loop form factors. The prefactors g?(A1l + 
BT) of each diagram consist of the identity and trace operator in flavor space 1 and T, 
respectively. For the BPS operator (1.2) and for the Konishi operator (1.4) the prefactors 
reduce to g*A and g?(A + BNg), respectively. They multiply the triangle integral which 
contains the numerator factor f (l). 


the numerator of the second integral in table 1 to squared momenta occurring in the 
denominator, it can be transformed to the expression 


(4 20— rn) =a) <<] --<X+ <L aK +h h-27) <<] l 


(H.5) 
Only the first three integrals are UV divergent. Moreover, they develop IR divergences 
if the corresponding external momentum square q?, p? or p3 vanishes. In this case, the 
respective integral vanishes identically in dimensional reduction since its IR pole and its UV 
pole cancel. The fourth integral is UV finite, but it becomes IR divergent if at least one of 
the three momentum squares vanishes. In case that p?, p3 are not zero, also the self-energy 
corrections of the scalar fields contribute to the form factor. Using the expression for the 
one-loop scalar self energy given in (H.4), the respective contribution can be written as 


KOKK K 


where the factor 5 originates from the fact that the square root of the renormalization 
constant determined from the self-energy contribution renormalizes the corresponding ele- 
mentary field. When added to the sum of the two diagrams given in table 1, this contri- 
bution exactly cancels the second and third term in the expansion of the second integral 
given in (H.5), irrespective of the vanishing or non-vanishing of p?, p2. In the case of the 
BPS operator, where both diagrams of table 1 only contribute with the coefficient A, the 
remaining UV divergence contained in the bubble integral cancels among the two diagrams 
given in table 1. Hence, in the BPS case, there is only a contribution from the triangle 
integral of (H.5). In the case of the Konishi operator, the contributions of the bubble 
integral do not cancel for the flavor-trace contribution, which comes with the coefficient B. 
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The one-loop form factors for the BPS operator and the Konishi operator hence read 


soso = (Bt + P3 — 2¢ = ' 
flo) = MCX + -2 =L 


We have calculated the above form factors for generic off-shell momenta p? 4 0 and p2 4 0. 


(H.7) 


Hence, they are generalizations of the respective expressions with p? = p2, given for the 
BPS operator in (3.4) and for the operator Ke in (3.16) into which the factor rg has to be 
introduced as prescribed in (4.4) in order to obtain the Konishi form factor. 

The contribution ft ) in (4.3) is free of any contribution from the triangle integral, 
and it is in particular independent of p? and pł. This explicitly confirms that the IR 
divergence is universal, i.e. the same for the BPS and the Konishi operator. Moreover, the 
UV divergence of the Konishi operator can be extracted from the final e-expansions of the 
Konishi and the BPS form factor given in (4.5) in the on-shell case p? = p2 = 0 where the 
the +-poles originate from both, the UV and the IR divergences. 


Two-loop form factors 


The one-particle-irreducible (1PI) diagrams for the two-loop form factors of the BPS oper- 
ator (1.2) and the Konishi operator (1.4) are displayed in table 2. The occurring integrals 
are given by 


E : dkar f(k,1) 

ta = Fd | = ome | Era 
7 i d?k dP”! f(k, 1) 

mte = (one | pe EEA 


E E kai f(k, D) 
hy = HK = ome | ST EENEN 
(H.8) 


For p? Æ 0, pł #0, contributions from diagrams involving the two-loop self-energy of 


the scalar fields have to be considered in addition to the 1PI diagrams shown in table 2 
Also, the second diagram coming with the integral J, yields a non-vanishing contribution, 
while it vanishes otherwise. All graphs are then IR finite, and the UV divergence can 
easily be extracted by setting e.g. one external momentum to zero and the other one to q? 
such that no new IR divergences are accidentally created. Moreover, since all integrals are 
superficially logarithmically divergent, one can neglect external momenta in the numerators 
as convenient for maximal simplifications. We have checked that this produces the known 
result for the two-loop overall UV divergence of the Konishi operator when subdivergences 
are subtracted by considering also the corresponding counter-term diagrams. This also 
produces a vanishing result for the BPS operator. 
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(k-I2(04 

(L+ p2)?(k +1) - (k +1 +29) 
(I+ q)? — 2k) - (T+ 2pa) 

(k —1)?(1 — p2) - (L+ 2q — p2) 


mai 
on 
x 
x 
x 
of (k — 1)? (l + 2p2) - (I+ 2p2 — q) 
x 
XK 
<Q 
< 
< 


(J + p2)? (l — 2k) - (1 — 2k — q) 


(1+ q)?(I+ p2)? 


PU a)? 


(k +1): (k +14 2q)(l — p2) - (l+ 2q — p2) 


- (L + 2p2) 


Table 2 — continued on next page 
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Iw 


Continued from previous page 
diagram g'(Al + BT) f(k, D 


(L+ p2)? - (k + 2q — p2) 


(l — 2k) - (l + 2p2)( 


(2k —1+ pa) - (l 

— (k +l + 2pz) - . Tot 

+ (21 — k + pr) - 

al - (—l — p2) 

— 2l- (k + 2q 

+2. (k—1I)(kṢ4 

I’ (k — pa): (k + 2q — pa) 


= rot + p2): (k — pe)(k + p2) - (k + 2q — p2) 


2l: (I+ p2)(k +q — p2) - k 
+21- k(l + po): (k +q -— p2) 


Table 2: Diagrams for the unrenormalized two-loop form factors. The prefactors g*(A1 + 
BT) of each diagram consist of the identity and the trace operator in flavor space, 1 and T, 
respectively. For the BPS operator (1.2) and for the Konishi operator (1.4), the prefactors 
reduce to gA and g*(A + BNg), respectively. They multiply the corresponding integral 
Ip, Iot or Ibp, which are given in (H.8) and contain the numerator factors f(k,1). For 
all diagrams which are not symmetric under a reflection at the horizontal axis, also the 
corresponding reflected version has to be considered. 


For p? = pa = 0, where the 1PI diagrams shown in table 2 are the only contributions 
to the form factors, it is advantageous to express the scalar products in the numerators in 
terms of squares of momenta found in the denominator from the propagators. Then, one 
can use IBP reduction as e.g. implemented in LiteRed [95] in order to further reduce the 
integrals to a set of master integrals. The results exactly match the ones given in (4.6). 
This confirms the absence of further rational terms that might not have been detected in 
the unitarity-based approach. 
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